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"Tenacity is one of the greatest skills or talents of a good engineer, particularly one who works in R&D. It is 
not unusual to work on-and-off on a problem for years before one feels like one has accomplished 

something noteworthy. 

The key to having such tenacity in the first place, though, is to believe in your own problem solving ability 
and to believe in the importance of the problem even in the face of doubt." 


Larry Young, NASA Ames Research Center 
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Abstract 

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, most of them 
are either unable to capture interference effects from bodies, or require an extremely large computational budget. 
The objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation 
of wind tunnel wall interference using a novel Computational Fluid Dynamics (CFD) solver for rotorcraft, RotCFD. 

In RotCFD, a mid-fidelity URANS solver is used with an incompressible flow model and a realizable k-c turbulence 
model. The rotor is, however, not modeled using a computationally expensive, unsteady viscous body-fitted grid, but 
is instead modeled using a blade element model with a momentum source approach. 

Various flight modes of the XV-15 isolated rotor, including hover, tilt and airplane mode, have been simulated and 
correlated to existing experimental and theoretical data. The rotor model is subsequently used for wind tunnel wall 
interference simulations in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames Research Center in 
California. 

The results from the validation of the isolated rotor performance showed good correlation with experimental and 
theoretical data. The results were on par with known theoretical analyses. In RotCFD the setup, grid generation and 
running of cases is faster than many CFD codes, which makes it a useful engineering tool. Performance predictions 
need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. 

For both test sections of the NFAC wall interference was examined by simulating the XV-1 5 rotor in the test section of 
the wind tunnel and with an identical grid but extended boundaries in free field. Both cases were also examined with 
an isolated rotor or with the rotor mounted on the modeled geometry of the Tiltrotor Test Rig (TTR). A 'quasi linear 
trim' was used to trim the thrust for the rotor to compare the power as a unique variable. Power differences between 
free field and wind tunnel cases were found from -7 % to 0 % in the 80- by 1 20-Foot Wind Tunnel test section and -1 .6 
% to 4.8 % in the 40- by 80-Foot Wind Tunnel, depending on the TTR orientation, tunnel velocity and blade setting. 
The TTR will be used in 201 6 to test the Bell 609 rotor in a similar fashion to the research in this report. 
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1 Introduction 

Experimental techniques to measure rotor and airframe aerodynamic performance are widely used but the need exists to 
understand the limitations of ground based testing by augmenting the analysis of experimental test results with 
Computational Fluid Dynamics (CFD) modeling. The objective of the present research is to develop an XV-15 Tilt Rotor 
Research Aircraft rotor model for investigation of wind tunnel wall interference. This research is performed to support wind 
tunnel tests scheduled for 2016. Ultimately the rotor model developed is used to investigate wind tunnel wall effects on 
large tilt rotors in the National Full-Scale Aerodynamics Complex (NFAC) facility at NASA Ames Research Center in California. 
The renewed interest in tilt rotors originates from NASA studies indicating significant reduction in congestion of 
commercial transport aviation 1 . 

The focus of this research is to understand the limitations and accuracy of tilt rotor performance predictions using a mid- 
fidelity CFD program. An unsteady RANS solver, RotCFD, is used with an incompressible flow and a k-£ turbulence model. 
RotCFD uses a model similar to an actuator-disk model (ADM) or blade-element model (BEM) with two-dimensional airfoil 
data that allows for relatively quick simulations of unsteady rotorcraft cases. The rotor is represented solely through the 
momentum it imparts to the flow. The coupling of the rotor with the surrounding flow is done by implementing its sources 
in the momentum equations. This omits the classical way of resolving a very fine grid around the rotor geometry to capture 
all flow effects. This gives RotCFD a significant advantage in simulation time 2 . 

Inherent to the use of two-dimensional airfoil data and a traditional blade-element model or actuator-disk model is the 
removal of three-dimensional effects. The most notable effects are stall delay, tip loss, yawed flow effects and unsteady 
rotor aerodynamics. Care must be taken these are properly accounted for either within RotCFD or by applying the proper 
models to the airfoil data tables in C81 format. 

The model will be applied to cases of the XV-15 rotor model on the Tiltrotor Test Rig (TTR) with struts in the two NFAC 
facilities under various flow conditions, pylon angles and rotor conditions. The performance with geometry, in terms of 
thrust and power, will be compared to isolated rotor performance in both free field and the wind tunnel test sections of the 
NFAC. 

7. 7 Research Question and Goal 

The main research question to be solved is phrased as follows: 

"To what extent will wind tunnel wall interference influence future performance measurements of the Civil 
Tilt Rotor (CTR) program in the National Full-Scale Aerodynamics Complex (NFAC) at NASA Ames Research 

Center?" 

The aim is to find the influence of the wind tunnel by modeling the TTR in the wind tunnel and in free field, and to 
determine the difference - if any - in performance. These deltas can be used to check the currently existing correction code 
and determine its compatibility with tilt rotor wind tunnel testing. The goal of this research is formally expressed as follows: 

"To understand and quantify the influence of wind tunnel testing to LCTR research programs at NASA 
Ames Research Center by assessing wind tunnel wall interference using a CFD approach. " 

The main goal is divided into sub goals that constitute the work breakdown structure as described in the project proposal 
and plan. 


1 W. W. Chung, D. Linse, A. Paris, D. Salvano, T. Trept, T. Wood, H. Gao, D. Miller, K. Wright, R. Young, and V. Cheng, "Modeling High-Speed Civil 
Tiltrotor Transports in the Next Generation Airspace," NASA CR 215960, pp. 1-414, 201 1 . 

2 Performance convergence for XV-1 5 isolated rotor in hover on baseline 2.5 Ghz. Intel i5 processor with 4 threads in less than half a day with 

implicit solver, steady model, 1 0 6 cells, and 1 0 3 * time steps. 
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2 Literature Study & Background 

Tilt rotors have the inherent capability to hover and fly in airplane mode or tilt mode. Because each of these modes has their 
own characteristics a short section of rotorcraft aerodynamics is included. This section will also provide a brief review on 
the Navier Stokes equations and CFD, and a short description of the National Full-Scale Aerodynamics Complex (NFAC) and 
Tiltrotor Test Rig (TTR). Lastly, the general characteristics of the XV-1 5 rotor are discussed and the selected experimental and 
theoretical performance data used for validation are presented. 

2.7 Rotorcraft Aerodynamics 

Hover is the defining characteristic of a helicopter; during hover the rotor produces lift but has no relative velocity to the air. 
A time-accurate representation of the wake and performance is very complicated to accurately compute. The analysis of a 
free field hover case also poses as the most computationally expensive as the rotor wake and inflow are solely developed by 
the rotor - no free stream component is present and thus relatively slow convergence occurs because of the slow induced 
velocity propagation. Furthermore, the boundaries must be set relatively far from the rotor in order to not influence the 
flow field. 

During edgewise forward flight the rotor moves through the air with a small forward tilt. This forward tilt contains the small 
forward component of the thrust vector that provides the propulsive force for the helicopter. In contrast to hover different 
flow phenomena arise here. Because of the freestream flow of air the advancing side of the rotor will experience a higher 
relative velocity than the retreating side of the rotor. This results in dissymmetry of lift and is normally accounted for by 
rotor trim. For this analysis no rotor trim will be considered and it is assumed that the total trust produced during trimmed 
forward flight is comparable to the untrimmed, 'averaged', rotor. The ratio of relative velocity over tip speed, V/QR, is used 
as definition for the advance ratio in this research. 

If the tilt rotor's pylon is lowered to airplane mode the rotor acts like an airplane propeller. Pylon angles between airplane 
and hover mode will be called tilt mode and can exhibit characteristics of both modes. Flight modes are considered with 
thrust values only well within the XV-15's flight envelope for this research to avoid (excessive) blade stalling. Autorotation 
will not be considered during this research. Furthermore, neither steady nor unsteady aeroelastic effects will be considered 
for this research as they are too small for the XV-1 5 to be important here [2]. 

2 . 7.7 Ground Effect Significance 

When a helicopter descends and reaches a proximity to the ground of around one rotor diameter it can experience ground 
effect. At the ground the wake velocity is reduced to zero which can result into pressure changes that can be transferred up 
to the rotor disk. With lower velocity and thus higher pressure on the bottom side of the disk the helicopter consumes less 
power for equivalent lift. This effect will increase as the distance to the ground is decreased. Figure 2.1 shows the general 
increase of thrust when the distance to the ground is decreased, according to various analytical models [3]. 



Ground effect influence decreases with forward speed due to deflection of the wake as can be seen in Figure 2.2. The 
ground effect is generally negligible when the rotor is more than one rotor diameter above ground [3]. 
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forward velocity over ideal hover induced velocity, V/v h M 

Figure 2.2 Ground effect; influence of forward speed on ground effect [3]. 

The graphs illustrating ground effect are only valid for use with a single ground plane. In a closed throat wind tunnel the 
complexity is increased as the rotor is effectively placed in a box with two open sides. 


2. 1.2 The Traditional Actuator Disk Model and Blade Element Model 


Several well-known models for rotors have been developed. The most well known are the Actuator Disk Model (ADM) based 
on momentum theory and the Blade Element Model (BEM). The traditional Actuator Disk Model is based on a uniform disk 
which supports and instantaneous pressure difference on either side of the disk. The pressure jump accelerates air over the 
disk until it reaches a finite velocity in the wake. Figure 2.3 shows the concept for an actuator disk in hover [4]. 



Flow Field Pressure Velocity 

Figure 2.3 Schematic representation of actuator disk in hover [4]. 


Assuming incompressible flow Bernoulli's equation can be evaluated before and after the rotor disk. Equation ( 2.1 ) shows 
the balance upstream and just before the rotor disk, while Equation ( 2.2 ) shows the balance just after the rotor disk and far 
downstream. 


Poo = Pi + 2 P V i 


( 2.1 ) 


1 1 

Pi+ Ap + -pv? =p m +-pv? (2.2) 

From Equation ( 2.1 ) and ( 2.2 ) the increase of pressure, Ap, is found to be 1/2 pv^. This change in pressure is, however, 
equal to the thrust per unit area, as shown in equation ( 2.3 ). Combined with the knowledge that, for momentum 
conservation, the thrust on the disc is the rate of increase of momentum, T = pAviv^ a second expression for Ap is obtained 
in equation ( 2.3 ). From Equations ( 2.1 ), ( 2.2 ) and ( 2.3 ) the relation for the downstream wake velocity according to 
momentum theory is now obtained in Equation ( 2.4 ). These values are used later on in the definition of the figure of merit. 


T pAv^ 

Ap = - = — ~ — = pvpv «, 


(2.3) 


v m = 2v { and T = 2pAv? 


(2.4) 
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Due to the homogenous disk the traditional actuator disk model is not very accurate when examining, for example, climb or 
edgewise flight. In these cases the non-uniform loading of the disk results in simulated performance values that can 
significantly diverge from the experimental ones. 

In the Blade Element Model (BEM) basic airfoil theory is applied to a rotating blade. Computation of the resultant velocity, 
taking into account the inflow, climb and rotational velocities, makes it possible to compute the sectional lift and drag 
forces. These are integrated to obtain the resulting thrust and torque coefficients. The accuracy was improved by 
implementing so-called prescribed or free wakes when the model was developed further. These wakes made it possible to 
generate the induced velocities at the rotor plane, yielding more accurate performance measurements. 

In RotCFD the model is placed in a CFD environment and the wake is developed by the NS equations. RotCFD incorporates 
two rotor models, both based on the sectional coefficients, which are placed in a CFD environment. The 'steady' model uses 
an averaged BEM to come to a model similar to the actuator disk. The term steady is chosen because within the disk the 
discrete blade loading is averaged. This is advantageous as this can yield less stringent time step restrictions. The disk is 
however not homogenous, in contrast to the 'traditional' ADM. The 'unsteady' model is based on a BEM, with discrete 
blades. 


From RotCFD the thrust, power and figure of merit are obtained internally. The thrust and power are non-dimensionalized 
using equations ( 2.5 ) and ( 2.6 ), respectively. The thrust is normal to the disk plane. 


T 

T ~ pA(Q.R) 2 


(2.5) 


The rotor shaft power is positive when power is supplied to the rotor. 

P 

Cp ~ pA(nR ) 3 


The geometric blade solidity is equal to the blade area over the rotor disk area. It is used for this calculation since only 
aerodynamic sections of the blade are modeled in this analysis (no hub or non-aerodynamic blade parts are included), as 
shown in Equation ( 2.1 0 ). The solidity value in literature is found to be 0.089 for the physical rotor [5]. The traditional value 
of this solidity, however, appears to assume that the chord extends all the way inboard to zero radius. The actual blade has 
both taper and root cutout [6]. 


NcR e _ NcR e 
A ~ nR 2 


0.081 


(2.7) 


The blade chord is assumed constant in this instance. The figure of merit is extracted directly from RotCFD. The full data 
table for the steady results is included in Appendix C. All plots for the unsteady model and full data table are included in 0. 


2.1.3 The Figure of Merit 

The figure of merit is a quantification of the efficiency of the rotor in hover. It compares the induced power, the power 
necessary for the lift obtained, to the total power, as shown in Equation ( 2.8 ). The total power consists of the induced 
power and the profile power, the power necessary to overcome the aerodynamic drag of the blades. 


M = 


Pi 

Pi+Po 




( 2 . 8 ) 


The last expression is found by substituting the profile power by non-dimensionalizing it with the rotor tip speed, DR, as 
shown in Equation ( 2.9 ). Furthermore the induced power is substituted with the use of Equation ( 2.10 ) and ( 2.1 1 ) to 
include the thrust coefficient. 


r = u 
p ° pA(D.R) 3 


(2.9) 
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The change in kinetic energy per unit time is 1/2 (pAv t )v^ which, combined with Equation ( 2.4 ) yields Tv t , or the 
definition of the induced power. Equation ( 2.1 1 ) shows how the thrust is non-dimensionalized using the representative tip 
speed. 


Pi = Tv t = T 3 ' 2 ^A 


( 2 . 10 ) 


T 

T “ pA(aR ) 2 


( 2.11 ) 


2. 1.4 Rotor Performance Corrections 

For a blade element or actuator-disk model there are several factors that might contribute to inaccuracies in the 
simulations. The most influential ones are: stall delay, tip loss, yawed flow effects and unsteady rotor aerodynamics that 
might all contribute to inaccuracies. RotCFD is assumed to internally resolve the last two effects due to the gridded rotor 
disk (or blade) and unsteady nature of the solver. Stall delay and tip loss must, however, be properly accounted for. The 
detailed approach to these corrections is described in Section 3.2. Because the C81 airfoil data tables and airfoil reference 
data is all acquired at full scale, no scaling corrections, or Reynolds number corrections, are adopted. The three 
dimensionality of the blade can, however, never be fully accounted for. 
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2.2 Navier-Stokes Equations and CFD 

CFD is the Analysis of fluid flow. One of the key drawbacks was the tremendous complexity of the behavior which let the 
economical use of CFD to be only possible for the wider industrial community since the 1990s [7]. With the rapid and 
continuous growth in computing power CFD research and development continues to improve efficiency and capability. 

The three important equations are the continuity, the momentum or Navier-Stokes equations, and the energy equations, 
which are in principal solely based on the conservation of mass, Newton's second law and the conservation of energy, 
respectively. RotUNS uses an incompressible flow, this means the energy equation can be decoupled form the analysis, and 
continuity and momentum equations are sufficient to solve for the velocity and pressure fields. The transfer of heat is not 
taken into consideration. 

2.2.1 Continuity Equation 

The continuity equation is based on the physical principle that mass can neither be created nor destroyed. The mass 
conservation is assured by setting the net mass flow out of the control volume through a surface equal to the time rate of 
decrease of mass inside the control volume. The continuity equation in conservation form is shown in Equation ( 2.1 2 ). 

^- + V-(pV) = 0 (2.12) 

Ot 

For an incompressible flow the density is constant so ^ = 0 and since the density is not equal to zero equation ( 2.1 2 ) can 
be rewritten to equation ( 2.1 3 ). 


V-V = 0 


( 2 . 13 ) 


2.2.2 Momentum Equation 

The momentum equations for a viscous flow are called the Navier-Stokes equations. The momentum equation is based on 
Newton's second law. The momentum equation in conservation form for the x-, y- and z-direction can be derived as shown 
in equation ( 2.14 ), ( 2.15 ) and ( 2.16 ), respectively [8]. Together they are called the Navier-Stokes equations, with 
representing the viscous stress tensor. 


d(pu) 

dt 


dp dr xx dr 


+ V ■ (pul 0 = - — + 

dx dx 


+ ■ 


yx 


dy 


9t zx 

dz 


+ Pfx 


( 2 . 14 ) 


d(pv) 

dt 


+ V 


f I r\ dp dr xy 



dj Z y 

dz 


+ Pfy 


( 2 . 15 ) 


d(pw) 

dt 


, N dp dx xz 
+5 . 0) “' l0= __ + _ 


d T V7 dTyy 

+ l7 + l? +p/z 


( 2 . 16 ) 


The left hand terms of these equations correspond to the unsteady and convection forces respectively. The terms on the 
right hand side correspond to the pressure, diffusion and external forces. 


By decomposing the NS equations (and the continuity equation) for incompressible fluids into a time-averaged part and a 
fluctuating part the Reynolds Averaged Navier-Stokes equations for incompressible flow can be derived. This approach still 
allows for unsteady simulations because only the random turbulence on small scale will be averaged out. The 
incompressible RANS equation is shown in Equation ( 2.1 7 ) [9]. 


dt 


+ 5j-(piW 


>--£♦ 


d(r 0 - - pu'tii'j) 

dxj 


( 2 . 17 ) 


2.2.3 Rotorcraft Computational Fluid Dynamics (RotCFD) 

The flow solver needs to handle unsteady rotor simulations while being computationally inexpensive. Rotorcraft CFD 
(RotCFD) is a mid-fidelity CFD tool specifically for rotorcraft design efforts and has been developed recently [10]. Young 
performed a study on complex rotor wake interaction simulation using RotCFD [11]. Rotorcraft Unstructured Solver 
(RotUNS) is a module within RotCFD that uses three-dimensional, unsteady, incompressible Reynolds-Averaged Navier- 
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Stokes equations (URANS) on a Cartesian unstructured grid with tetrahedral body-fitting near the body. The SIMPLER (Semi- 
Implicit Method for Pressure-Linked Equations Revised) is used in combination with the under relaxation factors to 
iteratively obtain the correct flow field in agreement with both the continuity and momentum equation. 

The module of RotCFD used in this research is RotUNS (Rotor Unstructured Solver). This name will be used synonymous to 
RotCFD from now on. Spalart discusses the possible pitfalls with the use of URANS but also indicates it is one of the few 
feasible unsteady methods when computational budget is limited [12]. laccarino et al. show the differences in simulated 
flow field using both a steady and unsteady Reynolds Averaged Navier-Stokes approach for the flow over a cube [1 3]. 

Turbulence is accounted for by the URANS equations combined with a two-equation realizable k-£ turbulence model with 
special wall treatment. The two transported variables are k, the turbulent kinetic energy, and £, the turbulent dissipation. 
Jones and Launder presented and validated the k-£ turbulence model in the 1970s [14], [15]. Yu & Cao have shown that 
accurate CFD analyses for the flow field and performance of a helicopter in forward flight can be obtained with a k-£ 
turbulence model with wall function method [16]. The wall function is introduced because the no-slip condition near the 
wall is found to behave unsatisfactory if the k-£ turbulence model is applied without this correction. 

Discretization is done using the finite-volume method and the implicit solver has been used in order to have a less stringent 
stability criterion and therefore more flexibility for various operations. This allows using consistent spatial and temporal grid 
under various cases more easily. The time-dependent solution, however, might not follow the exact transients as accurately 
as an explicit approach [1 7]. After the initial choice this has not been further investigated for this research. 

2.2.4 The RotCFD Rotor Model 

The rotor is modeled only through the momentum it imparts on the flow. "From the point of the fluid particles, the influence of 
the spinning rotor is to change their momentum." as stated in the original paper on this approach [18]. This momentum 
change occurs because of the aerodynamic forces exerted on the cells through the spinning blades. The fully viscous, 
unsteady, body-conforming grid usually chosen for such a problem is not necessary anymore, greatly reducing the 
complexity. Also the no-slip condition on the rotor is omitted and no unsteady boundary layer has to be resolved on a very 
dense grid, which also reduced the complexity. This reduction in complexity in turn is manifested through the possibility to 
run problems on desktop-class computers, where unsteady rotor simulations using full Navier-Stokes equations with 
turbulence models and body-conforming grid generally require supercomputers [18]. 

The source depends on the flow properties, rotor geometry, and two-dimensional airfoil data as shown in Equation ( 2.18 ). 
Where a and Cd are the airfoil section coefficients, a is the angle of attack, Vats is the absolute instantaneous velocity vector 
at the airfoil section, Q is the rotational velocity, x,y,z are the Cartesian coordinates, t is time, c is the chord, p is the density 
and b is the number of blades. 


5 = S(ci,c d ,a,V abs ,£l,x,y,z,t,c, p,b) (2.18) 

Figure 2.4 shows a schematic drawing of the rotor and the (r,0,z) inertial coordinate system fixed at the rotor used in the 
derivation of the momentum sources. 



Figure 2.4 Rotor and coordinate system [1 8]. 


The rotor model uses a version of the BEM called the 'steady' rotor model where the discrete rotor source terms can be 
averaged over the disc. This creates a model comparable to the traditional ADM, with the major differences that the disk is 
not, or doesn't have to be, experiencing a homogenous load. An 'unsteady' rotor model based on the blade-element model 
is also tested; this model does use the discrete blades, hence 'unsteady', blade modeling. The effect on the fluid is 
considered only at the point where the rotor, or disc if the 'steady' model is used, intersects with a grid point from the CFD 
domain. In the unsteady model a line at the quarter chord location replaces the chord length of the rotor. 
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Figure 2.5 shows the nomenclature of a blade section (see A-A in Figure 2.4) used to derive the source terms s</> and s z , 
aligned with tangential force coefficient, c f , and normal force coefficient, c n , respectively. The absolute velocity, Vats, consists 
of components v r , v</>and v z and the relative velocity, V re \, constists of components v/, v^'and v/. 



The normal and tangential (to the disk) coefficients can be derived from Figure 2.5 as shown in Equation ( 2.1 9 ) and ( 2.21 ). 

c n = c l cos/? — c d sin/? (2.19) 

c t = c x sin /? + c d cos / 3 ( 2 - 20 ) 

The angle /?, the angle between the relative Velocity, V re \, and the disk plane is obtained from Equation ( 2.21 ). Note that the 

relative velocity in 4> -direction is composed of the sum of the local flow component and the component through blade 

rotation at the radial station the blade section is chosen. 


v z 

tan /? = — (2.21) 

V’cp 

And the angle of attack, a , is now deduced in Equation ( 2.22 ). 

cc = 6-p (2.22) 

The section properties of the airfoil are determined from the angle of attack and Mach number, derived solely from the 
relative velocity, and inputted into Equation ( 2.19 ). The normal and tangential force coefficients can now be converted to 
4> and z-direction, as shown in equation ( 2.23 ) and ( 2.24 ). 

c z = -c n (2.23) 


The rotor blades are divided into spanwise elements. The blade geometric properties are constant over an element, and 
hence the source terms per unit element can be derived as shown in Equation ( 2.25 ) and ( 2.26 ). 


= 


S(p 


c z 


^„c) 


(2.25) 

(2.26) 
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These terms are changing as the rotor moves through the grid. In case the 'steady' model is chosen, the sources terms are 
averaged over the rotor disk while spinning. An example of the differences in the velocity field at the rotor disk for a steady 
and unsteady rotor model is shown in Figure 2.6 and Figure 2.7, respectively. Defining the time for one rotor revolution as 
trev and the time taken by a blade section to go through the width of a cell as t ce n and b the number of blades the averaged 
sources terms can be derived in Equations ( 2.27 ) and ( 2.28 ), assuming hover. For other flight modes the individual blades 
will have to be averaged as they experience different loading. 


bs z tcell ^^^z^cell 
tre v ^71 


( 2 . 27 ) 


bSq)tcell QbSqjtcell ( 2 . 28 ) 

trev ^71 



Figure 2.6 Velocity plot of a steady model. Figure 2.7 Velocity plot of an unsteady model. 


The source terms, defined in the inertial coordinate system fixed at the rotor, can be converted to the Cartesian coordinate 
system used in RotCFD. Integrating the z-coordinate of the source terms yields thrust whereas the integrated ^-direction 
source can be multiplied by the rate of rotation to obtain the (profile) power. The induced power is obtained by multiplying 
the thrust with the induced velocity. 

Please note the rotor precone angle has not been taken into account in this derivation, but is used in RotCFD and this 
research. However, no sweep angle is considered. Note how no component in r-direction is deduced, implying no radial 
effects are considered. Figure 2.8 and Figure 2.9 shows a representative example of the rotor (source) program grid with the 
cell structure underneath. The blades are only shown for visualization purposes. 



Figure 2.8 Top view of rotor disk grid and flow Figure 2.9 Enlarged view of the rotor disk grid and flow field grid, 
field grid. 
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2.2.5 Some considerations with the use ofRotCFD 

The RotCFD approach leads to an indirect relationship between the viscous effects and the action of the rotor blades. This 
approach therefore cannot give you the same results as the full NS solution would. Also dynamic stall, or any effect for 
which the boundary layers or actual rotor geometry must be known will probably not be evaluated properly. An example is 
the radial flow on the rotor blades that cannot be modeled. This is also the exact reason a stall delay model is evaluated in 
this research. 

Even though the rotor is based solely on two-dimensional airfoil data, because of the implementation into a flow field grid, 
some of the three-dimensional effects, for example tip loss, can be expected. RotCFD can incorporate flapping or cyclic 
pitch, but this is thought to be beyond the scope of the current research. 

2.3 National Full-Scale Aerodynamics Complex (NFAC) 

The isolated rotor wind tunnel tests for the tilt rotor research programs at NASA Ames Research Center are scheduled in 
May 2016 in the National Full-Scale Aerodynamics Complex (NFAC, Figure 2.10) 40- by 80-/80- by 120-Foot Wind Tunnels 
that are managed and operated by the U.S. Air Force. The original XV-15 research has been performed in the 40- by 80- 
Foot Wind Tunnel [1 9]. Zell has written elaborate reports on the performance and flow characteristics of both wind tunnels 
[20], [21]. Further aerodynamic characteristics are described by Corsiglia [22]. 



Figure 2.1 0 The NFAC facility at NASA Ames Research Center. 3 


2.3 . 7 Wind Tunnel and Test Section Geometry 

The two wind tunnels that comprise the NFAC facility share portions of their flow path. The cross section of the 80- by 1 20- 
Foot Wind Tunnel is an open circuit wind tunnel with a closed, rectangular test section (Figure 2.1 1). The 40- by 80-Foot 
Wind Tunnel is a single-return, closed section wind tunnel with an oval test section (Figure 2.1 1). 




Figure 2.1 1 Sketch of the 80- by 120-Foot Wind Tunnel and 40- by 80-Foot Wind Tunnel, respectively [20]. 


3 Anonymous, "The 80- by 1 20-foot wind tunnel, located at NASA Ames Research Center." 2005. [Online]. 
Available:http://www.nasa.gov/centers/ames/multimedia/images/2005/nfac.html. [Accessed: 23-Mar-2015]. 
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The geometry of the cross section is obtained from NFAC studies [20], [21] and sketched in Figure 2.12. The distinct 
differences in test section cross sectional geometry will result in a different flow behavior when the rotor flow is simulated. 
Because of the larger test section of the 80- by 120-Foot Wind Tunnel it is expected to show less interference with wind 
tunnel walls. The 40- by 80-Foot Wind Tunnel, besides having a smaller test section, is expected to show more wind tunnel 
interference with the rotor wake due to its smaller and oval cross sectional test section shape. The dimensions of the actual 
cross section are reduced slightly after the addition of acoustic liner. The inner dimensions of the test section are shown in 
Figure 2.12. The boundary layer profile was found not to change significantly with any of the tunnel operating variables 
[ 21 ]. 




Figure 2.12 The 80- by 1 20-Foot Wind Tunnel and 40- by 80-Foot Wind Tunnel cross section in the NFAC [20], [21]. 


2.3.2 Wind Tunnel Wall Corrections 

For positive thrust the velocity in the rotor slipstream will have a higher velocity than the wind tunnel velocity. Since both 
tunnels are closed-throat continuity dictates the velocity outside of the slipstream has a lower velocity. This lower velocity 
has a higher static pressure that works its way to the back of the rotor increasing the thrust [23]. The thrust could therefore 
result in a rotor reaching equal thrust at a lower free stream velocity. However, it might be possible under low-thrust (i.e. 
cruise) conditions, that the wake of the struts, with lower than free stream velocity, counteracts this effect in exactly the 
opposite way. Pope and Rae [24] describe the influence of rotor testing in tunnels based on earlier work by Glauert [23]. 
Figure 2.1 3 shows a schematic sketch of a rotor with wake in a wind tunnel and generic velocities. This sketch is only valid 
for airplane mode flight (axial flow). 


Velocity < V 



Figure 2.13 Propeller in closed throat wind tunnel [24]. 
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In edgewise flight this wind tunnel interference effect might be combined with an experienced 'ground effect' as the wake 
could reach the wall below a certain free stream velocity threshold. It is likely this effect is only to be expected in the 40- by 
80-Foot Wind Tunnel, if at all. During tilt mode, a combination of this effect with the axial flow interference effects could 
take place. All in all it is expected that pressure buildup will result in higher pressure differences on the rotor increasing 
performance in the wind tunnel, compared to true free field conditions. The characteristics of the NFAC which have been 
used in this research are presented in Appendix A. 



2.4 Tiltrotor Test Rig (TTR) 

During the actual wind tunnel tests in the NFAC, the isolated rotor will be mounted on the Tiltrotor Test Rig to allow for the 
control of powered rotor tests. The TTR (Figure 2.15 and Figure 2.16) is the modern version of the Proprotor Test Rig (PTR) 
after it needed replacement to allow for more advanced rotors. The TTR is mounted on three struts and four refurbished 
electric wind tunnel engines provide the power to the rotor. Table 2.1 contains the main design capabilities of the TTR. At 
the time of writing the calibration of the TTR has finished a couple of months ago and the stage of powered tests have been 
started. 



Figure 2.1 5 The TTR on the aft calibration stand, with cowlings open. 5 



Figure 2.1 6 The author at the TTR with the calibration rig mounted. 6 


4 Alexander Helmer, "The author at the 80-by 1 20-Foot Wind Tunnel Inlet.", 201 5. 

5 Anonymous, "TTR on the aft calibration stand, with cowlings open." 201 2. [Online]. 

Available: http://rotorcraft.arc.nasa.gov/Research/Facilities/ttrSlide.html. [Accessed: 24-Mar-201 5]. 

6 Alexander Helmer, "The author at the TTR with the calibration rig mounted", 2015. 
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Rotor forces will be measured on a dedicated balance installed within the TTR, instead of using the wind tunnel balance 
system. The wind tunnel turntable can rotate the TTR from axial to edgewise flight and all angles in between. The strut 
geometry in the 80- by 1 20-Foot Wind Tunnel has not yet been decided on at the time of writing [6] and will therefore be 
based on the 40- by 80-Foot Wind Tunnel struts. 

Table 2.1 Key TTR Design Capabilities [25]. 


Maximum Wind Speed (kts) 

Axial (Airplane) 

300 

Edgewise (Hover) 

180 

Rotational Speed (rpm) 

Minimal 

126 

Maximal 

630 

Maximum Thrust (lbs) 

Steady 

20,000 

Peak 

30,000 
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2.5 XV- 15 Rotor Characteristics 

The rotor parameters are summarized in Table 2.2. The value of the precone is that used for the Outdoor Aerodynamic 
Research Facility (OARF) test, which will be discussed later. The NACA 64-X25 airfoil at radial station r/R = 0.2500 [~] is copied 
to the root value at the cutout at r/R = 0.0875 H (marked with the in Table 2.2) to provide aerodynamic information at 
the root and prevent erroneous extrapolation of airfoil data. This, instead of using an extrapolation method, is common 
practice in rotor research [6]. The original airfoil sectional data files and adjusted or corrected files are public domain but 
have not been included because of their size. 

Table 2.2 XV-1 5 Rotor Characteristics [5], [6], [26]-[29]. 

Blade Geometry 


diameter 

7.62 m 

(25 ft.) 

disc area 

45.6 m 2 

(491 ft. 2 ) 

blade chord @.0875R 

0.432 m 

(17 in.) 

blade chord @1.000R 

0.356 m 

(Min.) 

blade area 

4.06 m 2 

(43.75 ft. 2 ) 

root cutout 

0.0875 r/R 


solidity 

0.089 


Blade Twist (bilinear) 

chord-line aerodynamic 

38.7° 


total chord 

41.5° 


Blade Airfoil Section, r/R [~] 

0.0875 

NACA 64-X25* 


0.2500 

NACA 64-X25 


0.5268 

NACA 64-X18 


0.8093 

NACA 64-X12 


1 .0000 

NACA 64-X08 


Rotor Characteristics 

hub precone angle 

2.5° 


Rotor rpm 

helicopter mode (hover, edgewise) 

589 


airplane mode (axial) 

517 


Blade tip speed 

helicopter mode (hover, edgewise) 

225.55 m/s 

(740 ft/s) 

airplane mode (axial) 

182.88 m/s 

(600 ft/s) 


2.5. 7 XV-1 5 Rotor Blade Characteristics 

The main parameters that are needed for the rotor are the twist and chord distribution, two-dimensional airfoil data along 
the span of the XV-1 5 blade, and the characteristic dimensions of the rotor. Some conflicting values were found in different 
publications [5], [26]— [29]; as an example the XV-1 5 was flown with multiple hub configurations leading to differently listed 
precone angles. The final values were decided in collaboration with experienced XV-1 5 researchers [6]. All input 
requirements for the computation of the rotor source, 5, in RotCFD are presented in Equation ( 2.29 ). 

S = S(c lf c d , a, V abs , H, x, y, z, t,c,p,b) ( 2.29 ) 
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The blade chord-line twist distribution is shown in Figure 2.1 7 and the blade chord distribution is shown in Figure 2.1 8. 



References to 'collective' will be implying the collective pitch angle, 0 O , which is defined zero when the pitch angle 
distribution equals the twist distribution in Figure 2.17. The relation with the quarter chord pitch angle, 0.75, is defined for 
the XV-15 rotor model in equation ( 2.30 ). 

0. 75 = 0 O + 6.61° (2.30) 



blade station, r/R [~] 

Figure 2.1 8 XV-1 5 rotor blade chord distribution [6]. 


Figure 2.19 shows the blade loci of the quarter chord, leading edge and trailing edge. The XV-15 rotor blade has 
approximately a -1 degree sweep angle that is not incorporated into the model because the sweep is only added for 
structural reasons [6]. 



blade station, r/R M 

Figure 2.1 9 XV-1 5 rotor blade aerodynamic loci [6]. 
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2.5.2 C81 Airfoil Data Structure 

The aerodynamic section properties for the blade airfoils elements are loaded from C81 data tables. The C81 airfoil data 
structure contains the airfoil section properties and originates from a first generation rotorcraft simulation program from 
Bell. The program is now outdated but the C81 format is still used, for example in CAMRAD II [30]. 

The format holds three two-dimensional arrays for the lift, drag and moment coefficients, Ct{(X,M), Cd{(X,M) and Cm{(X,M), 
respectively, as a function of angle of attack and Mach number. The format allows for an additional set for a trailing edge 
flap, this is however not used in this report or following analyses. The data format uses a separate rectangular array for each 
coefficient. The basic format is 10 columns, each 7 characters wide. The first column only holds the reference angle of 
attack, a, values. If the amount of Mach number entries is greater than 9, more than one line is used for each table row (a 
value). A new table row (a value) must start on a new line. An example of the partial data for the lift coefficient in C81 
format is shown in Figure 2.20. A brief summary of the format is further described in the CAMRAD II manual [30]. 


• • • ■ xvl 564x08_clean.c81 - 

(64-xea 2-0 AERODYNAMICS 122312231223 

0.0000 0.3000 0.4000 0.4200 0.5000 0.5200 0.6000 0.6600 0.7100 
0.7600 0.8100 0.8500 

-180,00-0. 1500-0. 1500-0.1500-0. 1500-0.1500-0. 1500-0.1500-0. 1500-0.1500 
-0.1500-0.1500-0.1500 

-90.00 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 
0.0000 0.0000 0.0000 

-45.00-1 . 0000-1 . 0000-1 . 0000-1 . 0000-1 . 0000-1 . 0000-1 . 0000-1 . 0000-1 . 0000 
- 1 . 0000 - 1 . 0000 - 1.0000 

-10.00-0. 7500-0.7400-0 . 6000-0.6000-0 . 5780-0. 5780-0 . 6120-0.7010-0 , 7000 
-0.6000-0.5500-0.4500 

-8.00-0 . 7400-0.7400-0 . 6000-0.6000-0 . 5780-0. 5780-0 . 6120-0. 7810-0 . 7200 
-0.7000-0.6000-0 . 5000 

-6.00-0.5410-0.5410-0.5540-0.5540-0.5780-0.5780-0.6120-0.7010-0.7600 

-0.7710-0.6800-0.5750 

-4.00-0.3400-0.3400-0.3530-0.3530-0.3680-0.3680-0.4010-0.4290-0.4820 

-0.5470-0.4400-0.4350 

-2.00-0. 1160-0. 1160-0.1200-0. 1200-0.1290-0. 1290-0.1450-0. 1540-0.1680 
-0.1930-0.1210-0.1150 

0.00 0.0910 0.0910 0.0950 0.0950 0.1010 0.1010 0.1110 0.1180 0.1300 
0.1500 0.2170 0.1850 


Figure 2.20 Example format of (partial) C81 airfoil data file. 

The C81 files are obtained from the Aeromechanics Branch at NASA Ames [6]. According to experienced XV-15 researchers 
the best available airfoil data set consists of four airfoils. The dataset is obtained using two-dimensional wind tunnel tests at 
full scale. The airfoil at radial station r/R = 0.2500 H is copied to the root value at the cutout, marked with the in Table 2.2. 


2.6 XV-15 Performance Data 

The XV-1 5 rotor will be used for this research because of the existing test data (wind tunnel and flight) and non-proprietary, 
publically available data. The XV-15 flight test data reports provide background in understanding the XV-15 and its 
performance [31]. 


XV-1 5 outdoor hover tests at NASA Ames have been documented by Felker and Betzina [28]. This data set was acquired on 
the Outdoor Aerodynamic Research Facility (OARF) test bed and will be referenced to as the 'OARF Data' further on in this 
report. This is believed to be the only full-scale data without wall effects for an XV-15 isolated rotor in hover. The Bell 
Helicopter company performed a study on the performance of an XV-15 isolated rotor in the NFAC facility [29] under 
various conditions, including edgewise and airplane mode. Johnson performed an assessment of the capability to calculate 
tilt rotor aircraft performance from this report [2]. It is unknown to what extent these results are influenced by the wind 
tunnel walls. 


Theoretical results obtained with CAMRAD I (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics) 
are also acquired from this reference. For the hover cases results from CAMRAD II [30], the improved version, are obtained as 
well. While sharing many characteristics with earlier versions of CAMRAD, CAMRAD II is completely recoded and has far 
more advanced options for computing rotor performance, loads, stability, etc. For predicting isolated proprotor 
performance, the most significant change is the introduction of stall-delay models. Here, the model developed by Corrigan 
and Schillings [32] is used. 

It is these two reports ([28] and [2]) and the hover performance from CAMRAD II [6], [30] that form the main basis for the 
validation of the rotor model to be developed. Because the V- 22 Osprey data is not publically accessible this is also the only 
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tilt rotor reference data available to validate the code. All performance data is obtained at sea level. The flight tests are 
assumed to be performed at a low enough altitude to render comparisons with data obtained from sealevel reasonable. 


2 . 6.7 Hover 

Figure 2.21 and Figure 2.22 show the power curve and figure of merit versus blade loading, respectively. The thrust and 
power are expressed as their non-dimensionalized coefficient values, divided by the rotor solidity. The theoretical power 
curves obtained from CAMRAD I and CAMRAD II show good agreement with various wind tunnel tests (WTT), flight data and 
the OARF data. Both curves and the corresponding data show an almost linear relation along the blade loading observed. 



blade loading, C T /o [~] 

Figure 2.21 XV-1 5 rotor hover power as a function of thrust [2], [6], [28]. 


The data points from the CAMRAD and OARF data are shown for the figure of merit plot in Figure 2.22. A considerable 
increase in scatter amongst the data points is observed. The figure of merit is defined as the ratio of induced over total 
power, as shown in Equation ( 2.31 ) [4]. 


M 


Pi 

Pi + Po 




1 + 


c Po VT 

r 3/2 

U T* 


(2.31 ) 


The peak of the curve consequently corresponds with a relative decrease of the induced power, and hence the total power 
is relatively increasing. This usually is the point where a significant section of the blade is stalling, and therefore the figure of 
merit, a measure for efficiency, drops. 



Figure 2.22 XV-1 5 rotor hover figure of Merit as a function of thrust [2], [6], [28]. 


It is thought the scatter along the figure of merit values is closely related to difficulties in obtaining clean power values in 
experiments or accurate drag values in theoretical analyses. 
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2.6.2 Tilt 

Several cases for a tilted rotor are observed in Figure 2.23. The pylon angle, a P , indicates the angle of the pylon, and rotor 
disk, to the relative wind velocity. a P = 0 [deg] corresponds to airplane or axial mode while a P - 90 degrees is helicopter or 
edgewise flight. The data is shown to deviate the most from theory at a pylon angle of a P = 75 degrees, attributed to 
deficiencies in the stall computation of CAMRAD I [2]. 


It is, however, crucial to note that the experimental values are obtained in the 40-by 80-Foot Wind Tunnel and might 
themselves exhibit wind tunnel interference effects. The tunnel interference was minimized by using various vents and 
panels that could be opened, reducing the blockage and alleviating the pressure. This approach has, however, not been 
further researched throughout this research. 


0.020 i 


0.016 


I 

o 

a 


0.008 


0.004 




/ 

/ 

/ 

/O 

/ 

✓ 

/ 

/ 

/ 

/ 



/ 

/ 

/O / 
/ / 

/ 

/o / 

/ 

/ 

/ 

/ 

/ 



/ , 
/ 

// 

/A 

\ 

\ 

\ 

\ 

\ 

\ 

V 

/ 

/ 

/' □ 

o 

/ 

o 

... 

/ A ,4 

, / 

/ 

□ ° 

° 


. 






0.000 


WTT ap = 1 5° (Mtip = .65) 

WTT ap = 60° (Mtip = .65) 

CAMRAD I ap = 15° (Mtip = .65) 

CAMRAD I ap = 60° (Mtip = .65) 


WTT ap = 30° (Mtip = .65) 

WTT ap = 75° (Mtip = .65) 

CAMRAD lap = 30° (Mtip = .65) 

CAMRAD I ap = 75° (Mtip = .65) 


0.00 


0.02 


0.08 


0.10 


0.04 0.06 

thrust, C T /a[~] 

Figure 2.23 XV- 1 5 rotor power as a function of thrust for different pylon angles at 1 //OR = .32 [2]. 

Figure 2.24 shows the influence of the advance ratio at a pylon angle of a P - 75 degrees. In this report the advance ratio is 
defined as the relative velocity magnitude over the blade tip speed, as shown in ( 2.32 ). 


V 

OR 


(2.32) 
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Figure 2.24 XV- 1 5 rotor power as function of thrust, for a p = 75° and Map = 0.65 [2]. 

The variation of advance ratio shows a consistent slope, but a diverging correlation as the advance ratio is increased. This 
deviation is thought to be attributable to the stall model in CAMRAD I, and can serve as a key comparison to rotor 
performance obtained with RotCFD. 


2.6.3 Airplane Mode 

Figure 2.25 shows the propulsive efficiency versus the blade loading for various tip speeds. The propulsive efficiency is 
calculated using Equation ( 2.33 ). Note the inherent difference to figure of merit, in Equation ( 2.31 ), which uses the 
induced velocity (evaluating useful power to obtain thrust) instead of the aircraft velocity (useful power to obtain 
propulsion). 




blade loading, C T /o [~] 

Figure 2.25 Rotor propulsive efficiency as function of thrust [2]. 


(2.33) 
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3 Aerodynamic Analysis Method 

This section describes the application and setup of the various validation models and wind tunnel models for this research. 

3.7 RotCFD Assumptions 

During this research RotCFD (RotUNS) is used to obtain rotor performance data, isolated or in a wind tunnel. The flow is 
considered incompressible, hence a Mach number in the whole flowfield must be lower than M = 0.3 for the flow to be 
considered low-subsonic or incompressible [33]. 

The rotor is solely modeled through the momentum it imparts on the flow. The rotor data is obtained from two-dimensional 
airfoil data. It is assumed this, combined with airfoil correction methods, will yield accurate rotor performance results. 
Because the rotor is modeled only through the momentum sources the tip Mach number is never truly present in the flow, 
only the momentum change in the surrounding cells. This makes the incompressibility assumption feasible. 

It is assumed that all Mach number and Reynolds number effects are considered through the momentum sources. Also 
dynamic stall effects, radial (boundary layer) flow effects cannot be properly obtained because there is no physical 
boundary layer, just the two-dimensional data. 

All simulations in this research are performed with SI units at sea level and International Standard Atmosphere (ISA). The 
values are summarized in Table 3.1. 


Symbol 

To 

Po 

po 

Oo 

R 

P 

Y 


Table 3.1 Sea level air properties used. 


Description 

Value 

sea level temperature 

288.16 [K] 

sea level air pressure 

101325 [Pa] 

sea level air density 

1 .225 [kg/m 3 ] 

speed of sound at sea level 

340.29 [m/s] 

gas constant 

287.05 [m 2 /s 2 K] 

dynamic viscosity 

1.75E-05 [m 2 /s] 

ratio of specific heats 

1.4 H 


3.2 XV- 1 5 Airfoil Data Corrections 

This section describes the airfoil data corrections performed on the XV-15 rotor model. For hover inclusion of unsteady 
aerodynamics and yawed flow aerodynamics has little effect on rotor thrust [32]. Furthermore it is assumed unsteady 
aerodynamics and yawed flow aerodynamics are covered to some extent in RotCFD's flow field computation. The other two 
main sources of inaccuracies are the stall delay effect and tip loss effect, both described in the following subsections. For 
forward flight unsteady aerodynamics and yawed flow effects are sufficient [32]. 

3.2.1 Stall Delay Model 

The effect of a rotor's rotation on the boundary layer is crucial for obtaining the correct performance prediction of a rotor 
[32], primarily during hover. Because the aerodynamic data files will provide the properties of the rotor, the boundary layer 
will not be resolved during the simulations. Therefore, the only way of taking into account the boundary layer rotational 
effect is altering the C81 airfoil data tables. 

Acree describes modeling requirements for analysis and optimization of the Joint-service Vertical take-off/landing 
Experimental Aircraft (JVX) proprotor performance [34]. A similar approach is used by altering the airfoil data tables using 
and applying stall delay according to the Corrigan and Schillings stall delay model [32]. The Corrigan stall delay model uses 
augmentation of the lift values in the C81 airfoil data tables by multiplying the section lift coefficient with a stall delay 
factor, K Li as shown in Equation ( 3.1 ). 


Radial flow along the blade span retards the point where the boundary layers breaks which delays stall, and therefore the 
maximum lift coefficient. Equation ( 3.1 ) shows the computation of the Corrigan stall parameter [35]. 

1 / 1 . 084 \ " 


K r = 


( C/T ( 

'.1517y 

v 136 V 

< c/r J 


= (1.291( c/r) 0775 )* 


(3.1 : 
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The exponent n varies from 0.8 to 1 .8, larger values usually giving better correlation [6]. Using n = 1 .8 for that reason the Kl 
values obtained are found in Figure 3.1 . 
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radial station, r/R [~] 

Figure 3.1 Spanwise Corrigan stall delay parameter obtained for XV-1 5 rotor model. 


Although the largest effect of stall delay is seen near the root, the tip also shows significant stall delay using this model. For 
Corrigan stall delay the stall delayed lift coefficient is now obtained in Equation ( 3.2 ). 


Cl = K l enable for CC < 30° 


(3.2) 


The stall delay is applied from 0 degrees angle of attack until 30 degrees angle of attack. From 30 degrees to 60 degrees the 
model is washed out using the washout, w, as displayed in Equation ( 3.3 ). 


w = 


60 — \a\ 
30 


(3.3) 


The lift for an angle of attack between 30 and 60 degrees can now be derived as displayed in ( 3.4 ). Over the angles of 
attack until 60 the lift coefficient approaches the regular a again. After 60 degrees no correction is applied anymore. 

Cl , washout = 0(tf L - 1) + 1) Q, table for 30° < a < 60° (3.4) 

The stall delay is only applied for positive section lift which is assumed to occur above an angle of attack of 0 degrees. It was 
found the zero lift angle of attack was close to, but not exactly equal to, zero and therefore small errors in the stall delay 
model might be introduced. These errors are however found to be negligible as the section lift within the angle of attack 
range from zero to zero-lift is very small. Comparison of the empirical model with experimental and theoretical data is 
presented in the original paper [32]. Section 4.5.2 discusses the results of the implementation of this model. 


3.2.2 Tip Loss Model 

Similar to aircraft wings, trailed vortex inflow over the tip of a rotor blade reduces its lifting capability. The RotCFD rotor 
model does not take this fully into account because it uses two-dimensional airfoil data and thus is likely to experience 
some section lift up to the blade tip. Flowever, the flow field environment of RotCFD will alter the relative velocities 
experienced by the blade. 


Leishman [36] describes an effective blade radius, or Prandtl tip loss factor, B, - usually around 98 % of the blade radius - 
that is unaffected by tip loss, as shown in Equation ( 3.5 ). This means the effective blade radius, R e , is a reduction of the 
actual blade radius, /?, (in terms of lift coefficient) based on the tip loss factor, B. The lift at the remaining 2 % of the blade is 
set to zero at the tip. Section 4.5.1 discusses the results of the implementation of the model. 

R e = RB (3.5) 


3.2.3 C81 Airfoil Adjustment Code 

The C81 data for the XV-1 5 airfoils vary almost linearly with radius, but the stall delay is non-linear, as shown in Figure 3.1. 
Because the stall delay is spanwise non-linear, the radial stations will have to be interpolated - even if the spanwise lift 
distribution is fairly linear. The airfoil data consists of the lift, drag and moment coefficients specified for a range of angles of 
attack and Mach numbers. A MATLAB code is written to perform a triple interpolation over angle of attack, Mach number 
and radial station. The interpolation over angle of attack and Mach number is necessary to make sure the radial 
interpolation can be properly evaluated. Once the angle of attack and Mach number matrix data is uniform for all three 
coefficients over the radial stations, the Corrigan stall delay model will be applied over a user-defined set of blade stations. 
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Figure 3.2 shows a representative set of coefficients versus angle of attack or Mach number imported for the tip airfoil, the 
NACA 64-X08 (C81 ) airfoil data file. For the other airfoils representative plots have been placed in 0. 





(d) c, vs. M at a = 4 [deg] for 64-X08 





Mach number, M [~] 

Figure 3.2 Angle of attack and Mach interpolation for representative cases of the NACA 64-X08 airfoil. 


For all airfoil files the angle of attack data is very well organized, showing high density of data points around the crucial 
areas. The Mach number data points are, however, less numerous. Especially for the root airfoils since they experience a 
lower relative Mach number. 


The green lines in the graph show the interpolation made by the code. Because of the numerous data points for the angle 
of attack the interpolation is done using a Piecewise Cubic Flermite Interpolating Polynomial (PCHIP) interpolation. This 
method was found to give the most 'natural' interpolation of rather unpredictable behavior of the curves. This method was 
not found to be working well for the interpolation of Mach number data points, here the PCHIP method could be 
unpredictable if one data point showed an abrupt deviation of the trend, as for example the case in Figure 3. 2d. Therefore 
the Mach number interpolation is done using a simpler linear interpolation. The following radial interpolation has also been 
done linearly. The interpolated lift curve is colored green in Figure 3.2. Figure 3.2a shows a clear stall region, albeit at a 
higher angle of attack than usually the case for aircraft wings, attributable to the relatively thick blades. This also manifests 
in the relatively gradual stall observed. The drop of the lift coefficient in Figure 3.2d shows the compressibility effects 
becoming troublesome. The spike in drag in Figure 3.2e and the drop in moment coefficient in Figure 3.2f show the results 
of the corresponding drag divergence, respectively. The stall delay and tip loss models only affect the lift coefficients and 
they will therefore be the focus variable in the remainder of the section. Figure 3.3 shows a simplified flowchart of the whole 
code for the airfoil adjustment, including the import, the interpolation along angle of attack, Mach number and radial 
stations, the stall delay and tip loss model implementation and the re-exporting to C81 format. 

The following steps are covered in the code: 

1. Start; 

2. Import of C81 airfoil data files and verification of the radial position of the airfoils; 

3. The C81 files are decomposed into matrices for easy manipulation; 

4. The angle of attack is interpolated using PCHIP-interpolation for all coefficients; 

5. The Mach number is linearly interpolated for all coefficients; 

6. The radial stations are interpolated for all coefficients; 

7. Application of stall delay, if selected in program; 

8. Application of tip loss factor, if selected in program; 

9. Conversion of matrices to C81 format and export to *.c81 files; 

10. End. 
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Figure 3.4 shows the spanwise lift coefficient of the original data, the data after angle of attack interpolation and Mach 
interpolation and the data after radial interpolation. The lift coefficient will mostly be discussed because both the stall delay 
and tip loss model only affect the lift coefficient. Note that the vectors specifying the remapping of the angle of attack, 
Mach number and radial stations are user specified. The plot is shown for a Mach number of M = 0.30 for continuous zero 
twist. Therefore these are not the actual lift coefficients the blade would exhibit because it is exhibits substantial twist as 
shown in Figure 2.17. 



Figure 3.4 Radial interpolation check of lift coefficient at M = 0.30. 

The results of the stall delay both for angle of attack and Mach number for various representative radial stations are shown 
in Figure 3.5. The effects of the model are clearly observed. Figure 3.5f is evaluated for a negative angle of attack, where no 
stall delay should be present. This is confirmed by the identical curves before and after model application. 


(a) c. vs. a at M = 0.50 [~] at r/R = 0.0875 [~] 




(c) c, vs. a at M = 0.71 [~] at r/R = 0.3000 [~] 



(d) c ; vs. Mat a = 30 [deg] at r/R = 0.0875 [~] 
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Figure 3.5 Stall delay influence as function of angle of attack or Mach number for various radial stations. 


More representative plots are presented in 0. 
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Tip loss is applied at r/R = 98% by setting the lift coefficient to zero after r/R = 98%. This drop occurs instantaneous [37] and 
is programmed to happen between r/R = .98 and r/R = .981. The final result before stall delay and tip loss compared to the 
unaltered interpolated data is shown in Figure 3.6. This plot shown for zero continuous twist at M = 0.30 and clearly shows 
the implication on lift coefficient of the application of both models. 



Figure 3.6 Lift coefficient as function of radial station at M = 0.30 at continuous zero pitch angle. 
The models will be tested independently during the validation of the rotor model in RotCFD. 


3.3 General Setup of Validation Cases 

The setup of the boundaries and gridding of the flow field is kept consistent through the validation cases. The boundaries 
consist of a rectangular prism with x,y,z-dimensions of 10D, 10D, 15D, respectively, with D being the rotor diameter. The 
rotor is placed in the center of the XY-plane and 5D in negative z-direction from the top plane. Note that positive z is 
aligned with thrust. This is done to eliminate any influence from the ground plane. The cell count is roughly between one 
and two million cells, depending on the case. The RotCFD user interface is shown in Figure 3.7 for a hover case; two side 
views of the applied gridding are shown in Figure 3.8 and Figure 3.9. 


N i i m O N W H 



Figure 3.7 RotCFD user interface with a hover case loaded. 
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3.3.1 Boundary Settings 

For a hover case all boundaries are pressure boundaries except for the plane in wake direction that is set to a mass outflow 
boundary correction. Cases where a free stream velocity is present have according velocity boundary corrections and a 
mass outflow boundary in the direction of the wake. 



Figure 3.8 Hover case, gridded centered side view (XZ) of the flowfield. 



Figure 3.9 Hover case, gridded top view (XY) at rotor height. 
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3.3.2 Spatial and Temporal Resolution Independency 

For all flight modes spatial and temporal resolution independency was checked by observing the performance parameters 
power, thrust and figure of merit to be independent of the chosen cell size and time step. A refinement box was used to 
confine the rotor in order to improve the result for the smallest cell count. The smallest cell size, found at the rotor disk 
itself, was found to be equivalent to roughly 8 [cm] whereas the time step for the validation with the steady model was set 
at 1/400 [s] and the time step for the unsteady model ranged between 1/800 [s] and 1/1200 [s]. For the unsteady model this 
corresponds to approximately 3 degrees of rotor rotation per time step. 

3.4 NFAC Wind Tunnel Cases Setup 

After the validation the interference in wind tunnels is examined for both wind tunnels of the National Full-Scale 
Aerodynamics Complex (NFAC) at NASA Ames Research Center. The two wind tunnels that comprise the NFAC facility share 
portions of their flow path. The cross section of the 80- by 120-Foot Wind Tunnel is an open circuit wind tunnel with a 
closed, rectangular test section. The 40- by 80-Foot Wind Tunnel is a closed section wind tunnel with an oval test section. 
For this research it is assumed the operating conditions are all fixed and at sealevel. 

The k-s turbulence model is primarily known for circulation in large areas, a wall function is added for more accurate 
velocity profiles near no-slip conditions at walls or other surfaces. No effort has been made to investigate the effect of 
chosen turbulence model versus a laminar approach or other turbulence models. 

Only the test section will be modeled with constant cross sectional area. In real life the test section of the open circuit 80- by 
120-Foot Wind Tunnel will always experience the influence from the turbulence from outside, the vane sets and the 
developed boundary layer and the flow therefore will therefore not be uniform. Similarly the test section of the closed 
section 40- by 80-Foot Wind Tunnel is likely to, despite the air exchange system, contain non-uniform flow because of 'old' 
disturbed air from the model. Disturbances can also occur at the test section inflow because of turbulent air from the tunnel 
walls, vane sets etc. To reduce the complexity of this research the inflow conditions are however chosen to be uniform at 
the test section inlet. 

The length of the test sections has been increased from their real value to eliminate influence on the rotor performance due 
to the imposed boundary conditions at the inlet and outlet. The test section elongation for both test sections in sketched in 
Figure 3.11, with the positive x-direction corresponding to the wind tunnel flow direction. For this research both test 
sections of the NFAC are therefore modeled as an open tunnel with closed test section. The pressure was found not to be 
disturbed around 2 to 3 rotor diameters distance from the inlet to the TTR. A render of the TTR in the 40- by 80-Foot Wind 
Tunnel test section during airplane mode and edgewise flight is shown in Figure 3.1 0. Note that this is the only direction the 
TTR can be rotated. 



7 E. Solis, "Render of TTR in 40- by 80-Foot Wind Tunnel Test Section." 2015. 
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Small features such as pressure tabs, small ramps or the faceting of the acoustic lining are assumed to have negligible 
influence on the results of this study and will not be modeled. The slight faceting on the 40- by 80-Foot Wind Tunnel will be 
ignored and the cross section is thus assumed to have semicircular walls. 

The XV-1 5 will be mounted on the Tiltrotor Test Rig (TTR). In 201 6 the TTR will be tested in the NFAC facility with the Bell 609 
rotor. Rotor forces will be measured on a dedicated balance installed within the TTR, instead of using the wind tunnel 
balance system; therefore the rotor performance from the rotor program in RotCFD can be used as the target variable. The 
wind tunnel turntable can rotate the TTR from axial (a P = 0 degrees) to edgewise flight (a P = 90 degrees) and all angles in 
between. Note that in edgewise flight the thrust will always be aligned with the positive y-axis and in airplane mode the 
thrust will always be aligned with the negative x-axis. 




Figure 3.1 1 The boundaries of the extended test sections with TTR in edgewise and axial mode, respectively. 


The thrust is manually trimmed using the collective setting to make sure values are similar to the validations results and the 
rotor is not evaluated far into the stall regime. 

RotUNS features a tetrahedral body-fitted grid with a Cartesian unstructured grid in the far field. RotUNS does not support 
viscous body-fitted grids and therefore will not accurately capture the boundary layer around geometry in the domain. The 
thickness of the boundary layers is found in Table A.1 and Table A.2 for the 40-by 80-Foot Wind Tunnel and the 80-by 120- 
Foot Wind Tunnel, but to reduce complexity it has not been tried to implement the effects of the actual tunnel wall 
boundary layer on the simulation. The choice for RotUNS was made because of the developmental stage of the RotVIS 
module at the time of writing and the experience that RotUNS should be mastered before attempting the use of RotVIS. 
Because the prime variable is rotor performance, not measured by forces on the body (i.e. TTR and/or struts), it is assumed 
the forces are of less importance to the predictions once performance convergence is guaranteed. 

A mass outflow condition is used at the end of the test section if positive wind tunnel speed is modeled. For hover cases a 
pressure boundary condition is applied to the inlet and outlet. The physical size of the grid relative to the rotor and in the 
region of the wake is kept identical to the validation cases. The time grid is kept equal to the one used in the validation 
cases, sometimes refined to keep stability at the tetrahedral body-fitted cells. Figure 3.12 shows two representative planes 
(XZ-plane or side view and YZ-plane or front view) of a fitted grid in both test sections. Note that the first (side) view is an 
extended test section to ensure uniform inflow due to the imposed boundary condition doesn't alter rotor performance. 
Lower grid density at the inlet is used to reduce cell count while the refinement at the walls is maintained, however, to keep 
an identical stream tube when the no slip condition at the walls results in velocity gradients over the adjacent cells. The 
total cell count for both test sections is between 900,000 and 1 ,200,000 cells. 
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80- by 1 20-Foot Wind Tunnel 


40- by 80-Foot Wind Tunnel 
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test section width: 24.38 [m] 


(extended) test section length: 62.94 [m] 

Figure 3.1 2 The extended test sections with TTR on struts and XV- 1 5 rotor in edgewise and axial mode, respectively. 



For each case evaluated there is a set of four tests consisting of the wind tunnel test section with geometry of the TTR and 
struts (WTGE) and without geometry of the TTR and struts (WTRO). Its free field counterpart without tunnel geometry but 
with (FFGE) and without the TTR and strut geometry (FFRO) is also modeled for direct comparison of rotor performance. The 
subsets for each case are summarized in Table 3.2. 


Table 3.2 Overview of four different subsets per case. 


Case 

Windtunnel 

Free Field 

Rotor 

TTR Geometry 

WTGE 

X 


X 

X 

WTRO 

X 


X 


FFGE 


X 

X 

X 

FFRO 


X 

X 



Rotor convergence is case dependent but occurs usually after the equivalent of around 10-15 rotor rotations and remains 
steady after that because the wake is out the influencing zone of rotor performance. This is observed visually by plotting the 
velocity and making sure the wake has progressed at least 1-2 rotor diameters from the rotor disk. Care must be taken that 
the rotor is effectively placed in a box (and not just in ground effect), which complicates the analysis. In some cases the 
wake might not have fully settled over the remainder of the geometry and thus force convergence on the TTR might not 
have occurred yet while rotor convergence is found. 

No automatic rotor trim is available yet in RotCFD at the time of writing. This means the thrust cannot be automatically 
maintained constant for each case. Instead the collective is fixed throughout the subsets. 

3.4 . 7 NFAC 80-by 120-Foot Wind Tunnel Cases 

The test section of the 80- by 120-Foot Wind Tunnel is expected to show little or no interference effect at all, or only under 
the most unfavorable conditions. Therefore axial flow at V = 1 00 [kts], the highest wind tunnel velocity, and edgewise at V = 
0 [kts] (hover) are the prime test cases. Hover will normally always be measured in axial mode for the least obstruction to 
the inflow and outflow of the rotor. However, the edgewise hover test is still evaluated to show the theoretical worst-case 
interference. An increase in tunnel velocity in edgewise mode should decrease the interference with the tunnel walls. 
Therefore the last case is an edgewise case at low tunnel velocity V = 10 [kts], expected to show smaller interference 
because of the higher tunnel velocity. Table 3.3 shows the three final cases for the 80-by 1 20-Foot Wind Tunnel test section. 

Table 3.3 Final 80- by 1 20-Foot Wind Tunnel Cases 


Case 

80- by 120-Foot Wind Tunnel 

V [kts] 

a P [deg] 

Oo [deg] 

Mti P H 

Airfoil Correction 

1 

edgewise 

0 

90 

4.00 

0.66 

TL, SD 

2 

edgewise 

10 

90 

4.00 

0.66 

TL, SD 

3 

axial 

100 

0 

14.80 

0.53 

TL 


Case 1, with an 'edgewise' hovering rotor, utilizes a pressure boundary condition at the in and outlet. Using the free field as 
a reference it is made sure the imposed pressure boundary is not altering the pressure changes altered by the rotor disk. 
Case 1 and 2 use collective settings equal to those used in the validation of the hover performance. Case 3 uses a collective 
setting roughly equivalent to low thrust, as would be expected for cruise, of around T ~ 3000 [N]. The WTGE cases for case 3, 
and case 1 and 2 are shown in Figure 3.13 and Figure 3.14, respectively. The wind tunnel velocity is aligned with the positive 
x-axis. 
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Figure 3.1 3 The setup of case 3, WTGE. Figure 3.1 4 The setup for case 1 and 2, WTGE. 


3.4.2 NFAC 40-by 80-Foot Wind Tunnel Cases 

The 40- by 80-Foot Wind Tunnel is likely to show interference under a larger set of conditions and will therefore be 
investigated in greater detail. Figure 3.15 shows the conversion corridor of the XV-15, which indicates the realistic operating 
conditions for the XV-1 5 rotor. This is checked to make sure no excessive stall or rotor mode simulation is asked for. 



For edgewise mode a similar case to the 80-by 1 20-Foot Wind Tunnel is investigated at V= 1 0 [kts] as well as a high velocity, 
V= 100 [kts], case. Figure 3.16 shows the height-velocity envelope of the XV-15. 
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It shows that the torque limit sets the cruise velocity at sea level at around V tas = 250 [kts]. The 40- by 80-Foot Wind Tunnel 
can run in excess of V = 300 [kts] but the incompressible solver will be limited to Mach = 0.3 [~] to avoid compressibility 
effects [33]. Therefore the axial mode at high speed is considered at V = 200 [kts]. The FFGE setup of cases 4 and 5 is 
pictured in Figure 3.17. The WTGE up for case 6 and the WTRO setup for case 7 is sketched in Figure 3.18 and Figure 3.19. 
Positive wind tunnel velocity is always aligned with the positive x-axis. 



Figure 3.1 7 The geometry for case 4 and 5, FFGE. 


Finally a tilted case with a P = 60 degrees is used to check the influence of the rotor disk advancing closer to the curved test 
section wall at V = 100 [kts]. The WTGE subsets of the edgewise and tilt mode are shown in Figure 3.18 and Figure 3.19, 
respectively. The final cases are summarized in Table 3.4. 
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Table 3.4 Final 40- by 80-Foot Wind Tunnel Cases 


Case 

80- by 1 20-Foot Wind Tunnel 

V [kts] 

o P [deg] 

Oo [deg] 

M tip [~] 

Airfoil Correction 

4 

edgewise 

10 

90 

4.00 

0.66 

TL 

5 

edgewise 

100 

90 

1.00 

0.66 

TL 

6 

axial 

200 

0 

32.68 

0.53 

TL 

7 

tilt 

100 

60 

8.00 

0.65 

TL 
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4 Validation of XV- 15 Rotor Model 

From RotCFD the thrust, power and figure of merit are obtained. The figure of merit is extracted directly from RotCFD. The 
full data table for the steady results is included in Appendix C. All plots for the unsteady model and full data table are 
included in 0. 

4.1 Hover 

The hover performance comparison for the unsteady model in terms of the power curve is shown in Figure 4.1. It shows 
serious deviation in excess of 10% from the reference data. The orange symbols show efforts in which several parameters 
were changed in an effort to find the cause of the bad correlation. In case El the grid count was increased from 1.3E6 to 
4.0E6 cells. In case E2 a spinner body was included, to check if the missing hub geometry showed a large influence because 
of root losses. Cases E3 and E4 had a significantly increased amount of iterations per time step and a lower relaxation factor, 
respectively. 
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Figure 4.1 Unsteady results for XV-15 rotor hover power as a function of thrust [2], [28]. 
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Figure 4.2 shows the same data points but expressed in terms of figure of merit. 


1.000 


0.800 


. 0.600 

u 

<U 

E 

o 

| 0.400 

O) 

u= 


0.200 


0.000 


0.00 





A < 
♦ 


A 

♦ 


cf> 

♦ 


o* 


A OARF Data (Mtip = 0.66) 

O WTT (Mtip = .53) 

CAMRAD I (Mtip = .66) 

RotCFD UNSTEADY TL (Mtip = .66) 

O RotCFD UNSTEADY TL SD E2 (Mtip = .66) 

O RotCFD UNSTEADY TL SD E4 (Mtip = .66) 

CAMRAD II (Mtip = .69) 


WTT (Mtip = .65) 

Flight test (Mtip = .69) 

RotCFD UNSTEADY (Mtip = .66) 

RotCFD UNSTEADY TL SD (Mtip = .66) 
RotCFD UNSTEADY TL SD E3 (Mtip = .66) 
RotCFD UNSTEADY SD (Mtip = .66) 


0.05 


0.10 

blade loading, C T /o H 


0.15 


0.20 


Figure 4.2 Unsteady results for XV-15 rotor hover figure of Merit as a function of thrust [2], [28]. 

The unsteady model again shows serious deviation from the reference data, generally around 10-20%. This error is very 
large and doesn't justify the use of the unsteady model. 


Figure 4.3 shows the results for XV-15 hover performance with the steady model. The correlation is generally good, 
although the tip loss (TL) model deviates from the reference data. The tip loss and stall delay model combined (TL SD) and 
the stall delay model (SD) both show very good correlation to both the OARF data and CAMRAD II curve. 
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Figure 4.3 Steady results for XV-15 rotor hover power as a function of thrust [2], [6], [28]. 
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Figure 4.4 shows the same data points for figure of merit. In this plot it is clear the stall delay only model has the best 
correlation with the OARF data points. Assuming the OARF data points are the most accurate dataset, the estimation of 
highest figure of merit is done most accurate using RotCFD. For the chosen tip speed not more OARF data points are 
available, but CAMRAD II shows serious stall effects as the curve steeply drops towards the higher thrust values. The stall 
delay only model clearly doesn't show a similar severe drop in performance towards the higher thrusts. 
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Figure 4.4 Steady results for XV-15 rotor hover figure of Merit as a function of thrust [2], [6], [28]. 

The steady model, in contrast to unsteady model, shows a fairly good correlation with the reference data and has therefore 
been chosen as the model for the wind tunnel simulations. The thrust and power as function of collective, for both the 
steady and unsteady model, are shown in Figure 4.5 and Figure 4.6, respectively, in an effort to find the source of the 
deviation of the unsteady model. It clearly shows the deviation in performance between both models, but no tested 
explanation has been found yet. 
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Figure 4.5 Rotor thrust as function of collective pitch angle. 
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Figure 4.6 Rotor power as function of collective pitch angle. 

Further work could include azimuthal check of lift, angle of attack and rotor conditions, to find out why the power is so 
much higher using the unsteady model. The lower figure of merit suggest a lower efficiency and thus more stall over the 
blade. 

4.2 Tilt Mode 

The evaluation of different pylon angles in tilt mode using the steady model is quite successful. No stall delay is applied on 
models except for hovering rotors. Figure 4.7 shows good correlation for all pylon angles up to (X P = 75 degrees. Noticeable 
is the difference with CAMRAD I results. The influence of tip loss factor is shown to be minimal. 
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Figure 4.7 Steady results for XV-1 5 rotor power as a function of thrust for various pylon angles at 1 //OR = .32 [2]. 
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Figure 4.8 shows the sensitivity to aspect ratio changes. The rotor model seems to perform well even up to the most severe 
aspect ratio tested. The influence of tip loss factor is shown to be minimal. 



blade loading, C T /o [~] 

Figure 4.8 Steady results for XV-1 5 rotor power as function of thrust, for a p = 75° [2]. 


Because the relatively low thrust and high aspect ratios the induced velocities at the tip are relatively smaller, therefore 
showing almost no influence of the tip loss factor, in strong contrast to the hover cases. Because the difference is so small 
the computational budget is chosen to be saved and not all cases are simulated for both models. 


4.3 Airplane Mode 


The final validation case is for two different tip speeds in axial, or airplane, mode. Correlation is very promising compared to 
the CAMRAD I curves. 
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Figure 4.9 Steady results for (rotor) propulsive efficiency as function of thrust [2]. 
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4.4 Ground Effect Study by Je Won Hong 

Je Won Hong from the University of Illinois at Urbana-Champaign offered his help during his internship at NASA Ames 
Research Center. He investigated, under guidance of the author, the XV-1 5 rotor, identical to the one used throughout this 
research, in ground effect. For this study he used the ideal power as this was the case for the theoretical reference models 
presented in literature [3]. Hong's work [38] shows that RotCFD shows expected behavior of rotors operating in ground 
effect as shown in Figure 4.1 0. 



4.5 Discussion on Airfoil Data Correction 

The original, tip loss model and stall delay model and their combinations have been examined in the present validation. The 
sections below include a discussion on the validity of the models and the obtained results. 

4.5.7 Tip Loss Factor 

It is likely that some of the three-dimensional effects of the tip (or root) of the rotor will be accounted for. Since the 
'conventional' tip loss model is applied, it would be expected the model overcompensates. As can be seen for both the 
steady and unsteady results the tip loss factor has a rather detrimental effect on figure of merit and power curve during 
hover. Especially the figure of merit plots in Figure 4.2 and Figure 4.4 clearly show this influence. The tip loss influence is 
considerably less during lower thrust modes, for example during airplane mode. 

There are a couple of factors that influence the tip loss; first of all the three-dimensionality of the blade tip aerodynamics or 
the induced velocity due to trailing vortices. Secondly because of the finiteness of the blade there is a pressure equalizing 
effect between the surfaces of the blade or rotor disk. 

The first mentioned reason is mostly present only in the unsteady model, whereas the second reason is applicable to both 
models. A further investigation of the tip loss model could yield better results. The main problem in RotCFD is thought to be 
the accurate computation of the angle of attack near the root or hub and its spanwise variation. Figure 4.1 1 shows a sketch 
for a typical observed difference in blade loading during hover between a traditional BEM, RotCFD and an actual blade 
loading. 


BEM 

RotCFD 

Actual 



radial station ► 

Figure 4.1 1 Sketch of blade loading characteristics. 
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4.5.2 Stall Delay Factor 

The stall delay factor is applied because of radial flow of the boundary layer that delays the onset of stall on the blade, 
particularly near the root of the blade. Since the rotor model described in Section 2.2.4 only accounts for normal and 
tangential source components, the radial direction is not taken into account. It would furthermore be very difficult, if not 
impossible, to create a similar system to also include radial boundary layers, because they would not be oriented in the 
direction of the main flow, and over relatively large lengths (the blade span) and therefore highly unpredictable. 

The Corrigan and Schillings stall delay model used is only valid for hover, for which is shows very good agreement with 
experimental and theoretical data. Although methods for (partial) stall delay in other flight modes exist, they are not 
pursued during this research. 

4.6 Accuracy and Precision 

An important factor in any simulation is the error observed. If there is a dataset that can be considered a benchmark the 
error can be calculated - this is, however, easier said than done. 

First of all, each data point represents a collective angle with a combined thrust and power variable. In past research, 
however, this collective angle has usually been neglected because to date no program is able to match collective setting 
with thrust and power accurately [6]. It can, however, produce an accurate data point on the power-thrust curve, despite a 
mismatch in collective angle. 

Without a way to directly compare data points, another possible way is to compare the RotCFD data points with a curve 
through the reference data. Nor RotCFD, nor the reference data, however, tend to have enough data points to justify a 
curve fit. The curve fit would also be subjective because no fitting technique is defined, or known. Another problem is that 
there is no way to know if the power, thrust or what combination of both generated the errors observed. 

The RotCFD data for the unsteady model shows good precision in general; any deviation from the reference data seems 
constant over the modeled thrust range. To get a first indication of the errors involved a very basic analysis has been 
performed. Figure 4.3 is repeated below in Figure 4.12, however, only including the RotCFD steady SD data points and the 
OARF data set. A second order polynomial fit has been added to the OARF data points. 

0.030 
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U! 0.020 
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J- 0.015 

| 0.010 

o 
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0.005 
0.000 

0.00 0.05 0.10 0.15 0.20 

blade loading, C T / o [-] 

Figure 4.1 2 2 nd order polynomial fit through the OARF data. 

By taking a representative thrust location on the curve and individually checking the error in power and the error in thrust, 
percentage errors can be obtained that serve as an indication of accuracy. These offsets where found to be around 3% for 
hover, 5% for tilt and 5% for airplane mode. The maximum figure of merit was 0.6% off from the value listed in the OARF 
data, which is exceptionally high, even for full NS codes, and is partially thought to be coincidentally this accurate. 

Of higher importance could be considered not the accuracy, but the precision. The precision is a measure of the consistency 
of the results and is associated with random errors. For the wind tunnel cases performance differences are obtained by 
finding the difference between two cases. This means the precision might be more important than the accuracy. The 
precision is, however, for the same reasons mentioned above, not obtainable without introducing new unknown errors into 
the data. The precision in the most important flight modes, hover, tilt and airplane mode, see Figure 4.3, Figure 4.7 and 
Figure 4.9, respectively is believed to be very good, as the data almost stays a consistent distance from the reference data. 
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4.7 Residual Values 

Most simulations for either hover, tilt or airplane mode show very similar characteristics in terms of residuals. Figure 4.13 
shows the moving average (MA) for the residuals of the main flow variables for a tilting isolated rotor at 15 degrees pylon 
angle, as shown in Figure 4.7. The length scale for the sampling of the moving average is the time-equivalent of one rotor 
rotation. The grey backgrounds show the original variation of the variables, while the colored lines show the moving 
average. This is done in order to not obfuscate all the overlapping values. The logarithm on the y-axis has base 1 0. 



It can be observed that momentum in all three directions shows the largest residual, or cumulative error over all the cells, of 
around 10°. From this number the location of highest errors in the flow field can unfortunately not be determined. The 
noticeable kink in the residuals slightly before time step 1000 is due to the wake of the rotor hitting the mass outflow 
boundary plane 1 0 rotor diameters downstream. Residuals can be seen to vary greatly with rotor wake development. 

Since the validation shows promising results it is concluded that the magnitude of the residuals, although relatively large 
sometimes, are not enough to offset the required accuracy for the rotor performance simulations. Increasing the iterations 
per time step or reducing the delta time will reduce the residuals, but no effect on rotor performance is observed. The 
residual overviews for hover, tilt at a 75 degree pylon angle and airplane mode are presented in Appendix C.2. 

4.8 Performance Convergence 

Performance convergence is found to be relatively robust. The thrust tends to converge rapidly, around the equivalent of 5 
rotor rotations, while the power tends to take a little longer. Usually around 10-15 rotor rotations the performance 
parameters, power and thrust for this research, have fully converged. To measure convergence a code has been written that 
plots the difference in performance compared to the final answer versus time steps; from this plot the percentage error can 
be determined which is used as indicator for convergence. If the difference is less than .25% to 1 % for a sustained amount 
of time, the solution is set to be converged. This amount of time is usually up to half of the simulated time. RotCFD has no 
internal convergence monitoring. 1 500 time steps roughly corresponds to 2 seconds of simulated time, or the equivalent of 
20 rotor rotations. 

Three performance convergence graphs are shown in Figure 4.14, Figure 4.15 and Figure 4.16, for hover, tilt and airplane 
mode, respectively. 
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Figure 4.1 4 Hover performance convergence over time. 
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Figure 4.1 6 Airplane mode convergence over time. 


Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics Witold J. F. Koning 


43 


Thrust, T [N] Thrust, T [N] Thrust, T [N] 


Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics Witold J. F. Koning 


44 


5 Results 

The wind tunnel cases required a substantial amount of time to run, varying from 2-6 weeks, depending on the required 
time difference between time steps for stability and convergence behavior of the case. Figure 5.1 shows a set of Apple 
iMacs used uninterrupted for the various subsets for about 2-3 months on end. 



Figure 5.1 Various computers used for months on end to compute each of the cases within the timeframe. 


5 . 1 NFAC 80-by 120-Foot Wind Tunnel Results 

The results for the rotor only cases and geometry cases for both free field and wind tunnel simulations are summarized in 
Table 5.1 and Table 5.2, respectively. The last column in both tables shows a top view (XY-plane) sketch of the wind tunnel 
configurations of the subset WTRO of each case. For the rotor only cases (FFRO, WTRO) the free field case (FFRO) is chosen 
as benchmark as this is the ideal testing ground. Similarly, the free field analogy with geometry (FFGE) is chosen as the 
benchmark in the geometry cases; FFGE and WTGE. Thus, the changes in power and thrust are expressed as relative 
percentage change compared to their free field counterpart. The computation of the change in power and change in thrust 
for FFRO and WTRO cases is shown in Equations ( 5.1 ) and ( 5.2 ). 

S p = - ~ PpFR0 xl00% (5.1) 


S t = 


t r FFR0 
— — — — — - xioo% 

^FFRO 


(5.2) 


Table 5.1 NFAC 80-by 1 20-Foot Wind Tunnel Rotor Only Results. 


Variable 

FFRO 

WTRO 

Sketch (WTRO) 

power, P [J/s] 

5.78E+05 

5.75E+05 

|t 

thrust, 7 [N] 

2.75E+04 

2.88E+04 

1 

change of power, <5 P [%] 

0.0 

-0.7 


change of thrust, <5 f [%] 

0.0 

4.7 

L> 

power, P [J/s] 

5.76E+05 

5.71 E+05 

V— t T 

thrust, 7 [N] 

2.77E+04 

2.89E+04 


change of power, <5 P [%] 

0.0 

-0.7 


change of thrust, <5 f [%] 

0.0 

4.3 

L* 

power, P [J/s] 

1.03E+05 

1.02E+05 


thrust, T [N] 

1.21 E+03 

1.19E+03 


change of power, <5 P [%] 

0.0 

-1.1 


change of thrust, <5 f [%] 

0.0 

-2.0 

L> 
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Case 1 and 2 show consistent reduction in power and increase in thrust, leading to a performance improvement of the wind 
tunnel cases compared to the free field cases. As expected the interference is reduced when a slight tunnel velocity of V = 
10 [kts] is used compared to edgewise hover. It is likely the interference will vanish after further increasing the tunnel 
velocity. Re-ingestion of the wake was studied for case 2 using a coarser model with an unfitted body. At V = 10 [kts] 
edgewise no re-ingestion of the wake occurred, as indicated in Figure 5.2. This 'coarse' run was simulated using a time step 
of around 1/40 [s] that is far below the validation temporal resolution. It does offer a rough idea of the wake propagation 
over longer time. 

Case 3 experiences performance change, but both power and thrust change in identical direction. This means the effect is 
likely to be less than observed for cases 1 and 2 because the data point can have moved along the thrust power curve. In an 
actual wind tunnel the rotor would have been trimmed for thrust. 


Table 5.2 NFAC 80-by 1 20-Foot Wind Tunnel Geometry Results. 


Case 

Variable 

FFGE 

WTGE 

Sketch (WTGE) 

1 

power, P [J/s] 

5.77E+05 

5.73E+05 

tl 


thrust, T[N] 

2.75E+04 

2.87E+04 

1 

f 


change of power, <5 P [%] 

0.0 

-0.7 

t 


change of thrust, <5 f [%] 

0.0 

4.6 

U. 

2 

power, P [J/s] 

5.75E+05 

5.70E+05 

V- . tr 


thrust, T[N] 

2.77E+04 

2.88E+04 

1 

j 


change of power, <5 P [%] 

0.0 

-0.9 

l 


change of thrust, 6 t [%] 

0.0 

4.0 

L* 

3 

power, P [J/s] 

1.39E+05 

1.38E+05 



thrust, T[N] 

1.95E+03 

1.93E+03 






T — ■ 


change of power, <5 P [%] 

0.0 

-0.5 



change of thrust, 6 t [%] 

0.0 

-1.2 

L. 


Differences between the cases in Table 5.1 and Table 5.2 are only apparent for case 3, in which the thrust and power of the 
case with TTR geometry is considerably higher compared to the rotor only cases, for an identical collective. While a true 
comparison of the two would need the thrust values to be trimmed, the difference, in part, is assumed to be attributable to 
an increase in pressure behind the rotor disc due to the inclusion of the wedge shaped front of the TTR. 



Figure 5.2 A velocity plot of a coarse, unfitted body test at t ~ 6 [s] showing no re-ingestion. 
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5 . 7.7 Quasi Trim for Thrust 

RotCFD lacks the functionality to trim the rotor automatically, at the time of writing. To further investigate case 1 -3 a 'quasi 
thrust trim' is performed for WTGE cases. By looking at the sign of the thrust change in the WTGE column, called WTGEi for 
now, a very small change in collective can be proposed to yield a second WTGE data point, called WTGE 2 . For this small 
change in P and T a linear interpolation is performed to find the power values of a WTGE data point with the same thrust as 
the FFGE data point. This 'quasi thrust-trimmed' WTGE data point is called WTGE*. This method allows showing the 
influence of the tunnel by isolating the power variable. A reduction in power suggests improved performance due to the 
wind tunnel, while an increase in power suggest a decrease in performance. The linearization is shown in Figure 5.3 and the 
computation is shown in Equation ( 5.3 ) and ( 5.4 ). 



Figure 5.3 Linearization to obtain point ( Tffge,P *). 


APwtpf P\a/tc Pwtffi 

^WTGE* — (TfFGE ~~ ?WTGEl) + ^WTGEl = CTt (TfFGE ~ ^WTGEl) + ^WTGEl ( 53 ) 

ni WTGE 'WTGE2 * WTGEI 

?WTGE* — ^FFGE ( 5A ^ 

The results of the linearization are summarized in Table 5.3. It can be seen the case 3 does not experience any noticeable 
tunnel interference. Case 1 and 2 both show and improvement in performance due to the tunnel geometry and as expected 
the performance difference decreases as tunnel velocity increases. The changes in power are, however, substantial, 
considering the size of the 80-by 120-Foot Wind Tunnel, primarily attributable to the influence of the walls as the 
differences with and without geometry are thought to be negligible. 


Table 5.3 Quasi Trimmed NFAC 80-by 1 20-Foot Wind Tunnel Geometry Results. 


Case 

Variable 

WTGEi 

wtge 2 

FFGE 

WTGE* 

1 

power, P [J/s] 

4.92E+05 

5.73E+05 

5.77E+05 

5.36E+05 


thrust, T[N] 

2.60E+04 

2.87E+04 

2.75E+04 

2.75E+04 


change of power, <5 P [%] 

- 

- 

0.0 

-7.2 


change of thrust, <5 f [%] 

- 

- 

0.0 

0.0 

2 

power, P [J/s] 

4.87E+05 

5.70E+05 

5.75E+05 

5.36E+05 


thrust, T[N] 

2.60E+04 

2.88E+04 

2.77E+04 

2.77E+04 


change of power, <5 P [%] 

- 

- 

0.0 

-6.7 


change of thrust, <5 f [%] 

- 

- 

0.0 

0.0 

3 

power, P [J/s] 

1.38E+05 

1.51 E+05 

1.39E+05 

1.39E+05 


thrust, T[N] 

1.93E+03 

2.19E+03 

1.95E+03 

1.95E+03 


change of power, <5 P [%] 

- 

- 

0.0 

0.0 


change of thrust, <5 f [%] 

- 

- 

0.0 

0.0 


For reference, the blockage fraction, expressed as the ratio of the frontal area of the TTR (without rotor) and struts to the 
wind tunnel cross sectional area, is found to be approximately 3.1% and 2.3% for the edgewise and axial case, respectively. 
These values are approximated using a CAD program. Representative velocity plots for case 2 are shown in Figure 5.4, more 
plots are presented in Appendix E. 
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p [Pa] : 9.69E+04 9.73E+04 9.76E+04 9.80E+04 9.84E+04 9.87E+04 9.91E+04 9.95E+04 9.98E+04 1.00E+05 1.01E+05 1.01E+05 1.01E+05 1.02E+05 1.02E+05 1.02E+05 1.03E+05 



V [m/s]:0.00E+00 5.11E+00 1.02E+01 1.53E+01 2.05E+01 2.56E+01 3.07E+01 3.58E+01 4.09E+01 4.60E+01 5.11E+01 5.62E+01 6.14E+01 6.65E+01 7.16E+01 7.67E+01 8.18E+01 

Figure 5.4 Representative velocity plots for case 3. 
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5.2 NFAC 40-by 80-Foot Wind Tunnel Results 

The results for the rotor only cases and geometry cases for both free field and wind tunnel simulations are summarized in 
Table 5.4 and Table 5.5, respectively. The last column in both tables shows a top view (XY-plane) sketch of the wind tunnel 
configurations of the subset WTRO of each case. 

Case 4 indicates low interference with both drops in thrust and power simultaneously. The proximity of the wall to the rotor 
in the 40- by 80-Foot Wind Tunnel seems to adversely influence the rotor at higher speeds as shown in case 5. A power 
increase with thrust loss is shown both with geometry and for the isolated rotor case. Care must be taken however, that the 
rotor in Case 5 is operating under a different collective and other free stream velocity, changing the blade loading and thus 
performance, making a comparison very difficult. 

Significant differences between isolated rotor and geometry cases are only observed in case 6. The difference is again 
thought to be caused by the increased pressure at the back of the rotor due to the wedge-like shape of the TTR nose, or by a 
difference in blade loading. This in turn increases the pressure in the vicinity of the TTR increasing the thrust. Case 6 with 
and without rotor shows a significant difference in performance. Due to constraints on the computational budget no 'quasi 
trim' has been performed on the rotor only cases, but it is thought that the operating conditions of the rotor are changed 
and likely shows a positive performance increase. Case 7 shows a thrust and power increase and therefore more favorable 
interference compared to the fully edgewise TTR in case 5. 


Table 5.4 NFAC 40-by 80-Foot Wind Tunnel Rotor Only Results 


Case 

Variable 

FFRO 

WTRO 

4 

power, P [J/s] 

4.84E+05 

4.79E+05 


thrust, T[N] 

2.47E+04 

2.44E+04 


change of power, <5 P [%] 

0.0 

-1.1 


change of thrust, 6 t [%] 

0.0 

-1.1 

5 

power, P [J/s] 

2.92E+05 

2.99E+05 


thrust, T[N] 

2.69E+04 

2.65E+04 


change of power, <5 P [%] 

0.0 

2.6 


change of thrust, 6 t [%] 

0.0 

-1.4 

6 

power, P [J/s] 

7.70E+05 

7.93E+05 


thrust, HIM] 

6.50E+03 

6.73E+03 


change of power, <5 P [%] 

0.0 

3.0 


change of thrust, 6 t [%] 

0.0 

3.6 

7 

power, P [J/s] 

7.26E+05 

7.44E+05 


thrust, HIM] 

2.28E+04 

2.37E+04 


change of power, <5 P [%] 

0.0 

2.5 


change of thrust, 6 t [%] 

0.0 

4.0 


Sketch 
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Similarly to the 80-by 120-Foot Wind Tunnel, re-ingestion has been examined for low tunnel velocities, i.e. case 4, and has 
been observed not to occur above at V = 1 0 [kts] edgewise. Similarly to the previous section, the geometry cases from Table 
5.5 are trimmed for thrust using a linearized method. 


Table 5.5 NFAC 40-by 80-Foot Wind Tunnel Geometry Results 


Case 

Variable 

FFGE 

WTGE 

4 

power, P [J/s] 

4.83E+05 

4.75E+05 


thrust, T[N] 

2.46E+04 

2.41 E+04 


change of power, <5 P [%] 

0.0 

-1.7 


change of thrust, 6 t [%] 

0.0 

-2.2 

5 

power, P [J/s] 

2.85E+05 

2.94E+05 


thrust, T[N] 

2.69E+04 

2.66E+04 


change of power, <5 P [%] 

0.0 

3.3 


change of thrust, 6 t [%] 

0.0 

-1.0 

6 

power, P [J/s] 

8.77E+05 

8.63E+05 


thrust, T[N] 

7.62E+03 

7.46E+03 


change of power, <5 P [%] 

0.0 

-1.6 


change of thrust, 6 t [%] 

0.0 

-2.0 

7 

power, P [J/s] 

7.30E+05 

7.46E+05 


thrust, T[N] 

2.33E+04 

2.41 E+04 


change of power, <5 P [%] 

0.0 

2.2 


change of thrust, 6 t [%] 

0.0 

3.1 


Sketch 



5.2.7 Quasi Trim 

The methodology to trim the thrust value is identical to the procedure in Section 5.1 .1 and Equation ( 5.3 ) and ( 5.4 ). The 
results of the quasi trim are show in Table 5.6. 


Case 

4 


Table 5.6 Quasi Trimmed NFAC 40-by 80-Foot Wind Tunnel Geometry Results 


Variable 

WTGE 1 

WTGE 2 

FFGE 

WTGE* 

power, P [J/s] 

4.75E+05 

5.08E+05 

4.83E+05 

4.88E+05 

thrust, HIM] 

2.41 E+04 

2.54E+04 

2.46E+04 

2.46E+04 

change of power, <5 P [%] 

- 

- 

0.0 

1.0 

change of thrust, 6 t [%] 

- 

- 

0.0 

0.0 

power, P [J/s] 

2.94E+05 

3.02E+05 

2.85E+05 

2.99E+05 

thrust, HIM] 

2.66E+04 

2.70E+04 

2.69E+04 

2.69E+04 

change of power, <5 P [%] 

- 

- 

0.0 

4.8 

change of thrust, 6 t [%] 

- 

- 

0.0 

0.0 

power, P [J/s] 

8.63E+05 

9.69E+05 

8.77E+05 

8.63E+05 

thrust, HIM] 

7.46E+03 

8.47E+04 

7.62E+03 

7.62E+03 

change of power, <5 P [%] 

- 

- 

0.0 

-1.6 

change of thrust, 6 t [%] 

- 

- 

0.0 

0.0 

power, P [J/s] 

7.46E+05 

4.96E+05 

7.30E+05 

7.25E+05 

thrust, HIM] 

2.41 E+04 

1 .50E+04 

2.33E+04 

2.33E+04 

change of power, <5 P [%] 

- 

- 

0.0 

-0.7 

change of thrust, <5 f [%] 

- 

- 

0.0 

0.0 


Both edgewise cases, case 4 and 5, show a decrease in performance, seemingly increasing with tunnel velocity. Both the 
axial and tilt case, case 6 and 7 respectively, show an increase in performance, as expected according to Glauert [23]. For 
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reference, the blockage fraction, is found to be approximately 7.7%, 3.8% and 7.5% for the edgewise, axial and tilt case, 
respectively. These values are approximated using a CAD program. Figure 5.5 shows the flow field of the FFGE subset for 
case 7, more flow field solutions are presented in Appendix E. 



Figure 5.5 Flowfield of FFGE subset for case 7. 


5.3 Performance Convergence 

The performance convergence for case 2 and 6 is shown in Figure 5.6 and Figure 5.7. Performance convergence below .25% 
is usually quickly observed, following the same procedure as for the validation cases. Case 2 has a high amount of time 
steps and shows Typical' slower convergence for hover cases. Case 7 shows rapid convergence characteristic for airplane or 
tilt mode cases. 
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3.50E+04 
3.00E+04 
2.50E+04 
2.00E+04 
1 .50E+04 
1 .00E+04 
5.00E+03 
O.OOE+OO 


5.4 Forces on the TTR 

The rotor performance will be directly measured within the TTR axis. Therefore the wind tunnel scale doesn't have to be 
used and the forces on the TTR are, for this research, not of much importance as RotCFD outputs the rotor performance 
directly. Moreover, since RotUNS doesn't allow for viscous body-fitted grids, chances of an accurate drag measurement 
greatly diminish as no plausible boundary layer will be modeled, or better, sustained. Also the unsteady nature of the 
URANS solver may show slight variations in the forces due to unsteady (vortex) shedding. 

Two representative graphs with the forces on the TTR are outputted, however, to show particular problem of interest. 
Figure 5.8 shows the integrated forces on the TTR for an edgewise case at high tunnel velocity - case 5, WTGE. The behavior 
of these forces over time is as expected, flattening out as the flow and wake are fully established over time. 

2.50E+04 


2.00E+04 

1 .50E+04 

Z 

LL. 

oT 1.00E+04 

u 
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o 

LL. 

5.00E+03 
0.00E+00 
-5.00E+03 

0.0 1,000.0 2,000.0 3,000.0 4,000.0 5,000.0 6,000.0 7,000.0 

timestep 

Figure 5.8 Integrated forces on the TTR for case 5, WTGE. 

A particular case, which has been run for extra time to investigate its behavior, is that of case 4. The edgewise low velocity 
case in the 40-by 80-Foot Wind Tunnel test section. It shows that even after nearly doubling the time steps the forces on the 
TTR have not converged yet. 
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Figure 5.9 Integrated forces on the TTR and struts for case 4, WTGE. 


The power and thrust only vary .5 % and 1% respectively over the last 50% of simulated time, as is shown in Figure 5.10. This 
extremely stable performance behavior of the rotor was found consistently throughout this report. It recommended for 
future work to further investigate the transient nature of the forces observed in this case. 
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Figure 5.1 0 Performance convergence for case 4, WTGE. 
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5.5 Stability Related Time Step Restrictions 

A problem encountered during this research was the stability related time step restriction. The 'clean' rectilinear Cartesian 
grid, such as that one used in the isolated rotor cases for the validation, was combined with the tetrahedral cells for the 
body fitting. The tetrahedrals show higher aspect ratios per cell and sharp angles, compared to the near-squared cells in the 
far field. In the case of stability related time step restrictions it is thought the turbulent kinetic energy, k, and turbulent 
dissipation values, e, diverge over time, which, as they are coupled back into the RANS equations, cause flow field 
divergence. 

The addition of these cells therefore limit the maximum time step used, often much smaller than the time step necessary for 
accurate rotor performance values as used in Section 0. In some cases a time step equivalent to half a degree of rotor 
rotation would be necessary, tremendously increasing the computational effort. 

It is left for future work to investigate the effect of another gridding technique, perhaps using RotVIS or another module and 
investigate the performance results running with more refined unfitted bodies solely using a Cartesian grid. Figure 5.11, 
Figure 5.12 and Figure 5.13 show the TTR nose ungridded (only meshed), with the used body fitting, and an unfitted 
version, respectively. 



Figure 5.1 1 Meshed TTR nose (no grid). 


Figure 5.1 2 Body fitted TTR nose. 


Figure 5.13 Gridded TTR nose without body 
fitting. 
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6 Conclusions and Recommendations 


6.1 Conclusions 

The setup and grid generation in RotCFD is faster than many CFD codes and therefore makes it a useful engineering tool. 
The validation results show very good correlation with performance data over the thrust range examined. The results for 
the XV-15 isolated rotor in hover with a Corrigan stall delay method showed excellent performance correlation compared 
the OARF data set. The tip loss model was found to overcorrect, most likely caused by the generated tip vortex or pressure 
equalization over the blade tip. The tip loss is, however, not fully accounted for by RotCFD, as can be seen from the 
obtained blade loading. For the tilt and airplane mode cases the XV-15 rotor model showed accurate results compared to 
wind tunnel data, with only minor differences observed between the tip loss model and clean variant. These differences are 
due to the lower thrust settings and lower induced velocities at the tips mostly because of the advance ratio. Overall the 
author believes a very good compromise between accuracy and efficiency is achieved. For hover tests the Corrigan and 
Schillings stall delay model with n = 1.8 is advised and for tilt and airplane mode the clean XV-15 rotor will suffice for very 
reasonable performance predictions in complicated flows. 

For the 80-by 1 20-Foot Wind Tunnel no influence was found for the XV-1 5 rotor in axial mode at the highest tunnel velocity 
of V = 100 [kts]. For the case with (near) hover in edgewise mode, the performance shows a power decrease of -7.2% and - 
6.7% for hover and edgewise flight with V = 10 [kts], respectively. The large interference at this very unfavorable 
configuration is attributed to the boxed-in effect of the rotor facing the tunnel wall. 

For the 40-by 80-Foot Wind Tunnel interference was found under all tested conditions as is to be expected due to the 
smaller cross sectional area of the test section. In edgewise mode both at tunnel velocities of V = 10 [kts] and V = 1 00 [kts] 
an increase in power is observed of 1.0% and 4.8%, respectively. Differences between the two edgewise cases are hard to 
justify because of the different operating conditions of the rotor. Both the axial case at V= 200 [kts] and the tilted case at V 
= 100 [kts] shows a decrease in power of -1.6% and 0.7%, respectively. This is expected for a thrusting rotor in airplane mode 
in a duct. The axial case at V = 200 [kts] shows a wake velocity slightly over M = 0.30, which could be questioned for its 
validity due to compressibility effects. Researchers at Ames were not convinced it would yield any substantial difference in 
this case, especially since the delta is obtained between two similar cases. 

6.2 Recommendations 

The investigation of the effect of stall delay in other flight modes is left for future work. A better tip loss model, either by 
reducing the tip loss factor or finding a better suitable method is also left for future work. At the time of writing a new 
model has been introduced for use with RotCFD. The approach corrects the local angle of attack to achieve zero lift at the 
tip based on the local pressure difference between the top and bottom of the rotor. 

All performance validation data is expressed as function of thrust coefficient over solidity. A better picture of the correlation 
is obtained if this solidity value is not used, and hence the coefficients only are used as parameters. This leads to more direct 
comparisons of the thrust and power. Because of the ambiguous solidity data (not every report is clear on what value is 
used) the solidity value is incorporated in the present calculations. 

Considerable amount of work can be done to investigate the blade loading of the XV-15 rotor. Because of the changing 
blade loading at different azimuthal stations and numerous flight modes and thrust values this is a considerable effort. The 
thrust modeled in this research is well in the normal range of the XV-15. More research can lead to performance data into 
the rotor stall region. The assumption that the rotor was not trimmed, to reduce the rolling moment, might also cause 
adverse rotor effects. Therefore cyclic pitch variation and its effect on performance can also be further investigated. 

RotUNS is capable of adaptive grids. For this research it was chosen to keep the gridding and flow field cell sizes nearly 
identical over the cases for consistent performance values and comparison. A tailored adaptive grid, however, could 
improve the wake propagation through the domain and due to the specified high cell density could yield a more efficient 
simulation or a more accurate result in equal computation time. It is highly interesting if an adaptive grid can yield 
improved computed velocities at the blade or disk tips and hence discard the need for a tip model at all. 

The reason for the performance mismatch using the unsteady model is not found. One hypothesis is that, even at a 
considerably small delta time, the implicit method fails to follow the exact time variations of the flow field. When this 
mismatch is then incorrectly matched with the 'hard coded' timing of the discrete rotor blades the time of the blade 
position and wake state are not correct and could lead to different pressures over the rotor disk. Also the rotor trim could 
affect the unsteady results more aggressively compared to the steady model. 
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A substantial issue found with the use of RotUNS and body fitted grids was the stability related time step restriction. 
Because of the need to reduce the time step considerably the computational time per case increased significantly. A finer 
grid doesn't necessarily solve the problem because the effective cell shape at the body fit remains similar. An investigation 
into the use of RotVIS, which can handle body-fitted viscous grids, and therefore should simulate boundary layers more 
accurately as well, might present useful as the body fitted cells are much more organized. The use of RotVIS for this research 
was highly discouraged before the author was fully familiar with RotUNS. If the computation time can be considerably 
reduced using RotVIS it might be possible to generate enough performance data points per case to actually generate power 
curves. These are much more useful to compare than individual data points. Another area of interest is the wall pressure 
along the longitudinal axis of the tunnel test sections. These can be obtained by extracting the pressure along the 
centerlines of the walls, presumable excluding the floor-plane. A second set of simulations for empty test sections must be 
performed, however. 
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Appendix A NFAC Characteristics 

The main characteristics of the NFAC wind tunnel test sections are summarized in Table A.1 and Table A.2. 

Table A.1 40x80 Foot Wind Tunnel Characteristics [21]. 


40 x 80 Foot Wind Tunnel 



Width Test Section 

24.38 m 

(80 ft) 

Height Test Section 

12.19 m 

(40 ft ) 

Length Test Section 

24.38 m 

(80 ft) 

Actual Width Test Section 

24.08 m 

(79 ft) 

Actual Height Test Section 

11.89 m 

(39 ft) 

Actual Length Test Section 

24.38 m 

(80 ft) 

Approximate BL Thickness Floor (Start - End) 

.25 - .46 m 

(10-18 in) 

Approximate BL Thickness Top (Start - End) 

.08 -.15 m 

(3 - 6 in) 

Approximate BL Thickness Sides (Start - End) 

> .25 - .46 m 

(> 10 -18 in) 

Maximal Test Section Velocity 

154.33 m/s 

(300 kts) 

Table A.2 80 x 1 20 Foot Wind Tunnel Characteristics [20]. 


80 x 120 Foot Wind Tunnel 

Width Test Section 

36.58 m 

(120 ft ) 

Height Test Section 

24.38 m 

(80 ft) 

Length Test Section 

36.58 m 

(120 ft ) 

Actual Width Test Section 

35.97 m 

(118 ft) 

Actual Height Test Section 

23.93 m 

(78.5 ft) 

Actual Length Test Section 

36.58 m 

(120 ft ) 

Approximate BL Thickness Floor (Start - End) 

.76 -1.12 m 

(30 -44 in) 

Maximal Test Section Velocity 

51.44 m/s 

(100 kts) 
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Appendix BC81 Airfoil Adjustment Code Results 


B. 7 Angle of Attack and Mach number interpolation 

The C81 airfoil files are obtained from experienced XV-15 researchers [6]. The interpolation for angle of attack and Mach 
number is performed. Figure B.1 shows part of the imported data for the 64-X08 airfoil and the PCHIP interpolation for 
angles of attack and linear interpolation for Mach numbers. 





(d) c, vs. M at a = 4 [deg] for 64-X08 





Figure B.1 Imported C81 data for NACA 64-X08 airfoil. 


Similar data is shown in Figure B.2, Figure B.3 and Figure B.4 for the 64-X1 2, 64-X1 8 and 64-X25 airfoils, respectively. 
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(e) c d vs. M at a = 2 [deg] for 64-X1 2 
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Figure B.2 Imported C81 data for NACA 64-X1 2 airfoil. 
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(b) c d vs. a at M = 0.60 [~] for 64-X1 8 




(d) c. vs. M at a = 4 [deg] for 64-X1 8 



(e) c d vs. M at a = 2 [deg] for 64-X1 8 
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Figure B.3 Imported C81 data for NACA 64-X1 8 airfoil. 
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(a) c vs. a at M = 0.75 [~] for 64-X25 







B.2 Representative Stall Delay Plots 

Figure B.5 and Figure B.6 show the effect of stall delay on the lift curve slope on a set of interpolated airfoil data at various 
radial stations. 






(e) c. vs. a at M = 0.85 [~] at r/R = 0.9700 [~] 



Figure B.5 Effect of stall delay on lift curve slope. 


(f ) c. vs. a at M = 0.85 [~] at r/R = 1 .0000 H 
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Figure B.6 Effect of stall delay on lift coefficient for set of airfoil data at various radial stations. 


Figure B.6d shows the result of the assumption that the zero lift angle of attack is equal to zero. As a result the lift 
coefficient is slightly lowered at high Mach numbers. While the observed effect is small, future work could eliminate this 
error from the program. 
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Appendix C Steady XV-i 5 Rotor Validation Results 

All the plots for the steady results are presented in Section 0. 

C. 7 Simulation Parameters 

Table C.3 shows a compact overview of the simulation parameters. ND stands for non-dimensionalized, the boundary size is 
expressed as the coordinates of the corners of a rectangular prism. The grid cells indicate the amount of cells on the x,y and 
z edge of the boundary, respectively. Rotor Grid Ref. and Grid Refinement indicate refinement of cells by multiplying the 
amount of cells by a certain factor. 


Table C.3 Overview of steady simulation parameters. 


Hover Parameters 

Boundary Conditions & Grid 

Simulation Variables 


Radius 

3.81 m 

Boundary Size 

(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.1) 

Simulated time 

2.5 s 

Cone angle 

2.5 deg 

Grid cells [#] 

30,30,42 

Timesteps 

3000 

Tip speed 

225.55 m/s 

Rotor Box Ref. 

1,1 

AT 

2.50E-03 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 

6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement 
Box [m] 

(-7.62, -7.62, -68.58) 
(7.62, 7.62, 38.1) 

Iterations 

10 

Power ND 

6.41E+08 

Grid Refinement 

4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

2.84E+06 

Volume Ratio Max 

8 

Flight Condition 

hover 



Cells 

1,053,948 

Rotor Model 

steady 

Tilt and Advance Ratio 

Boundary Conditions & Grid 

Simulation Variables 


Radius 

3.81 m 

Boundary Size 

(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.10) 

Simulated time 

1.25 s 

Cone angle 

2.5 deg 

Grid cells [#] 

30,30,42 

Timesteps 

1000 

Tip speed 

221.17 m/s 

Rotor Box Ref. 

1,1 

AT 

2.50E-03 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 

6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement 
Box [m] 

(-15.24, -7.62, -68,58) 
(7.62, 7.62, 38.10) 

Iterations 

10 

Power ND 

6.05E+08 

Grid Refinement 

4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

2.73E+06 

Volume Ratio Max 

8 

Flight Condition 

general 



Cells 

2,029,692 (at 30 degree tilt) 

Rotor Model 

steady 

Airplane Mode 


Boundary Conditions & Grid 

Simulation Variables 


Radius 

3.81 m 

Boundary Size 

(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.10) 

Simulated time 

1.25 s 

Cone angle 

2.5 deg 

Grid cells [#] 

30,30,42 

Timesteps 

1000 

Tip speed 

183.76 m/s 

Rotor Box Ref. 

1,1 

AT 

2.50E-03 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 

6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement 
Box [m] 

(-15.24, -7.62, -68,58) 
(7.62, 7.62, 38.10) 

Iterations 

10 

Power ND 

3.47E+08 

Grid Refinement 

4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

1.89E+06 

Volume Ratio Max 

8 

Flight Condition 

general 



Cells 

1,053,948 

Rotor Model 

steady 
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C.2 Residual Overview 

The present section shows the residuals for a representative airplane mode, tilt mode and hover mode case in Figure C.7, 
Figure C.8 and Figure C.9, respectively. 


2 r 



500 1000 

timestep 

Figure C.7 Residual overview for representative airplane mode case. 


1500 



timestep 

Figure C.8 Residual overview for representative tilt mode case with ap = 75. 
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6 r 



1500 

timestep 

Figure C.9 Residual overview for representative hover mode case. 


C.3 Data Steady XV- 7 5 Rotor Results 

The data used in the various plots is tabulated in Table C.4. 

Table C.4 Summary of steady RotCFD Validation Data 


# 

Flight Mode 

File 

a P 

[deg] 

V/OR 

[m/s] 

00 

[deg] 

TL 

SD 

T[N] 

P [J/s] 

Ct la [-] 

Cp/G [-] 

M [-] or 

nH 

1 

Hover 

O.rpr 

90 

0.00 

0.0 

N 

N 

1.467E+04 

2.487E+05 

6.373E-02 

4.790E-03 

0.675 

2 

Hover 

2.rpr 

90 

0.00 

2.0 

N 

N 

1.952E+04 

3.483E+05 

8.480E-02 

6.708E-03 

0.741 

3 

Hover 

4.rpr 

90 

0.00 

4.0 

N 

N 

2.457E+04 

4.799E+05 

1.067E-01 

9.243E-03 

0.759 

4 

Hover 

6.rpr 

90 

0.00 

6.0 

N 

N 

2.972E+04 

6.474E+05 

1.291E-01 

1.247E-02 

0.738 

5 

Hover 

8.rpr 

90 

0.00 

8.0 

N 

N 

3.466E+04 

8.578E+05 

1.506E-01 

1.652E-02 

0.710 

6 

Hover 

lO.rpr 

90 

0.00 

10.0 

N 

N 

3.867E+04 

1.122E+06 

1.680E-01 

2.161E-02 

0.640 

7 

Hover 

0_TL.rpr 

90 

0.00 

0.0 

Y 

N 

1.450E+04 

2.485E+05 

6.299E-02 

4.786E-03 

0.665 

8 

Hover 

2_TL.rpr 

90 

0.00 

2.0 

Y 

N 

1.909E+04 

3.464E+05 

8.293E-02 

6.672E-03 

0.718 

9 

Hover 

4_TL.rpr 

90 

0.00 

4.0 

Y 

N 

2.373E+04 

4.810E+05 

1.031E-01 

9.264E-03 

0.720 

10 

Hover 

6_TL.rpr 

90 

0.00 

6.0 

Y 

N 

2.844E+04 

6.612E+05 

1.235E-01 

1.273E-02 

0.685 

11 

Hover 

8_TL.rpr 

90 

0.00 

8.0 

Y 

N 

3.282E+04 

8.967E+05 

1.426E-01 

1.727E-02 

0.678 

12 

Hover 

10_TL.rpr 

90 

0.00 

10.0 

Y 

N 

3.643E+04 

1.180E+06 

1.583E-01 

2.273E-02 

0.557 

13 

Hover 

0_TL_SD.rpr 

90 

0.00 

0.0 

Y 

Y 

1.669E+04 

2.906E+05 

7.250E-02 

5.597E-03 

0.703 

14 

Hover 

2_TL_SD.rpr 

90 

0.00 

2.0 

Y 

Y 

2.183E+04 

4.060E+05 

9.483E-02 

7.819E-03 

0.752 

15 

Hover 

4_TL_SD.rpr 

90 

0.00 

4.0 

Y 

Y 

2.712E+04 

5.617E+05 

1.178E-01 

1.082E-02 

0.753 

16 

Hover 

6_TL_SD.rpr 

90 

0.00 

6.0 

Y 

Y 

3.245E+04 

7.686E+05 

1.410E-01 

1.480E-02 

0.721 

17 

Hover 

8_TL_SD.rpr 

90 

0.00 

8.0 

Y 

Y 

3.759E+04 

1.031E+06 

1.633E-01 

1.986E-02 

0.669 

18 

Hover 

10_TL_SD.rpr 

90 

0.00 

10.0 

Y 

Y 

4.207E+04 

1.338E+06 

1.828E-01 

2.577E-02 

0.609 
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19 

Tilt 

15_18_TL.rpr 

15 

0.32 

18.0 

Y 

N 

8.109E+03 

6.146E+05 

3.663E-02 

1.255E-02 

NA 

20 

Tilt 

15_19_TL.rpr 

15 

0.32 

19.0 

Y 

N 

1.208E+04 

9.008E+05 

5.457E-02 

1.840E-02 

NA 

21 

Tilt 

15_20_TL.rpr 

15 

0.32 

20.0 

Y 

N 

1.600E+04 

1.193E+06 

7.227E-02 

2.436E-02 

NA 

22 

Tilt 

15_21_TL.rpr 

15 

0.32 

21.0 

Y 

N 

1.985E+04 

1.496E+06 

8.966E-02 

3.055E-02 

NA 

71 

Tilt 

15_175_TL.rpr 

15 

0.32 

17.5 

Y 

N 

6.133E+03 

4.765E+05 

2.770E-02 

9.731E-03 

NA 

72 

Tilt 

15_185_TL.rpr 

15 

0.32 

18.5 

Y 

N 

1.009E+04 

7.559E+05 

4.558E-02 

1.544E-02 

NA 

23 

Tilt 

30_155_TL.rpr 

30 

0.32 

15.5 

Y 

N 

8.998E+03 

5.399E+05 

4.064E-02 

1.103E-02 

NA 

24 

Tilt 

30_165_TL.rpr 

30 

0.32 

16.5 

Y 

N 

1.289E+04 

7.906E+05 

5.822E-02 

1.615E-02 

NA 

25 

Tilt 

30_175_TL.rpr 

30 

0.32 

17.5 

Y 

N 

1.673E+04 

1.046E+06 

7.557E-02 

2.136E-02 

NA 

26 

Tilt 

30_185_TL.rpr 

30 

0.32 

18.5 

Y 

N 

2.045E+00 

1.312E+06 

9.237E-06 

2.679E-02 

NA 

73 

Tilt 

30_15_TL.rpr 

30 

0.32 

15.0 

Y 

N 

7.030E+03 

4.174E+05 

3.175E-02 

8.524E-03 

NA 

74 

Tilt 

30_16_TL.rpr 

30 

0.32 

16.0 

Y 

N 

1.095E+04 

6.644E+05 

4.946E-02 

1.357E-02 

NA 

27 

Tilt 

60_8_TL.rpr 

60 

0.32 

8.0 

Y 

N 

1.512E+04 

5.360E+05 

6.830E-02 

1.095E-02 

NA 

28 

Tilt 

60_9_TL.rpr 

60 

0.32 

9.0 

Y 

N 

1.865E+04 

6.911E+05 

8.424E-02 

1.411E-02 

NA 

29 

Tilt 

60_10_TL.rpr 

60 

0.32 

10.0 

Y 

N 

2.207E+04 

8.546E+05 

9.969E-02 

1.745E-02 

NA 

30 

Tilt 

60_ll_TL.rpr 

60 

0.32 

11.0 

Y 

N 

2.525E+04 

1.027E+06 

1.141E-01 

2.097E-02 

NA 

75 

Tilt 

60_6_TL.rpr 

60 

0.32 

6.0 

Y 

N 

7.823E+03 

2.571E+05 

3.534E-02 

5.250E-03 

NA 

76 

Tilt 

60_7_TL.rpr 

60 

0.32 

7.0 

Y 

N 

1.149E+04 

3.904E+05 

5.190E-02 

7.972E-03 

NA 

31 

Tilt 

75_2_TL.rpr 

75 

0.32 

2.0 

Y 

N 

1.374E+04 

3.304E+05 

6.206E-02 

6.747E-03 

NA 

32 

Tilt 

75_3_TL.rpr 

75 

0.32 

3.0 

Y 

N 

1.716E+04 

4.160E+05 

7.751E-02 

8.495E-03 

NA 

33 

Tilt 

75_4_TL.rpr 

75 

0.32 

4.0 

Y 

N 

2.041E+04 

5.117E+05 

9.219E-02 

1.045E-02 

NA 

34 

Tilt 

75_5_TL.rpr 

75 

0.32 

5.0 

Y 

N 

2.358E+04 

6.180E+05 

1.065E-01 

1.262E-02 

NA 

77 

Tilt 

75_l_TL.rpr 

75 

0.32 

1.0 

Y 

N 

1.030E+04 

2.575E+05 

4.652E-02 

5.258E-03 

NA 

35 

Advance Ratio 

75_2_TL.rpr 

75 

0.32 

2.0 

Y 

N 

1.374E+04 

3.304E+05 

6.206E-02 

6.747E-03 

NA 

36 

Advance Ratio 

75_3_TL.rpr 

75 

0.32 

3.0 

Y 

N 

1.716E+04 

4.160E+05 

7.751E-02 

8.495E-03 

NA 

37 

Advance Ratio 

75_4_TL.rpr 

75 

0.32 

4.0 

Y 

N 

2.041E+04 

5.117E+05 

9.219E-02 

1.045E-02 

NA 

38 

Advance Ratio 

75_5_TL.rpr 

75 

0.32 

5.0 

Y 

N 

2.358E+04 

6.180E+05 

1.065E-01 

1.262E-02 

NA 

81 

Advance Ratio 

75_l_TL.rpr 

75 

0.32 

1.0 

Y 

N 

1.030E+04 

2.575E+05 

4.652E-02 

5.258E-03 

NA 

39 

Advance Ratio 

27_75_2_TL.rpr 

75 

0.27 

2.0 

Y 

N 

1.534E+04 

3.345E+05 

6.929E-02 

6.831E-03 

NA 

40 

Advance Ratio 

27_75_3_TL.rpr 

75 

0.27 

3.0 

Y 

N 

1.857E+04 

4.099E+05 

8.388E-02 

8.371E-03 

NA 

41 

Advance Ratio 

27_75_4_TL.rpr 

75 

0.27 

4.0 

Y 

N 

2.165E+04 

4.965E+05 

9.779E-02 

1.014E-02 

NA 

82 

Advance Ratio 

27_75_l_TL.rpr 

75 

0.27 

1.0 

Y 

N 

1.209E+04 

2.673E+05 

5.463E-02 

5.459E-03 

NA 

42 

Advance Ratio 

18_75_3_TL.rpr 

75 

0.18 

3.0 

Y 

N 

2.053E+04 

3.982E+05 

9.273E-02 

8.132E-03 

NA 

43 

Advance Ratio 

18_75_4_TL.rpr 

75 

0.18 

4.0 

Y 

N 

2.325E+04 

4.704E+05 

1.050E-01 

9.606E-03 

NA 

44 

Advance Ratio 

18_75_5_TL.rpr 

75 

0.18 

5.0 

Y 

N 

2.596E+04 

5.519E+05 

1.173E-01 

1.127E-02 

NA 

86 

Advance Ratio 

18_75_l_TL.rpr 

75 

0.18 

1.0 

Y 

N 

1.490E+04 

2.792E+05 

6.730E-02 

5.702E-03 

NA 

87 

Advance Ratio 

18_75_2_TL.rpr 

75 

0.18 

2.0 

Y 

N 

1.771E+04 

3.349E+05 

8.000E-02 

6.839E-03 

NA 

45 

Airplane Mode 

40_23_TL.rpr 

0 

0.40 

23.0 

Y 

N 

2.906E+03 

2.762E+05 

1.902E-02 

9.837E-03 

0.777 

46 

Airplane Mode 

40_24_TL.rpr 

0 

0.40 

24.0 

Y 

N 

5.753E+03 

4.902E+05 

3.765E-02 

1.746E-02 

0.866 

47 

Airplane Mode 

40_25_TL.rpr 

0 

0.40 

25.0 

Y 

N 

8.577E+03 

7.101E+05 

5.613E-02 

2.529E-02 

0.892 

48 

Airplane Mode 

40_26_TL.rpr 

0 

0.40 

26.0 

Y 

N 

1.140E+04 

9.343E+05 

7.461E-02 

3.328E-02 

0.901 

78 

Airplane Mode 

40_225_TL.rpr 

0 

0.40 

22.5 

Y 

N 

1.480E+03 

1.706E+05 

9.686E-03 

6.076E-03 

0.640 

79 

Airplane Mode 

40_255_TL.rpr 

0 

0.40 

25.5 

Y 

N 

9.988E+03 

8.215E+05 

6.537E-02 

2.926E-02 

0.898 

49 

Airplane Mode 

70_375_TL.rpr 

0 

0.70 

37.5 

Y 

N 

3.328E+03 

5.535E+05 

2.178E-02 

1.971E-02 

0.777 
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50 

Airplane Mode 

70_385_TL.rpr 

0 

0.70 

38.5 

Y 

N 

7.193E+03 

1.054E+06 

4.707E-02 

3.754E-02 

0.882 

51 

Airplane Mode 

70_395_TL.rpr 

0 

0.70 

39.5 

Y 

N 

1.098E+04 

1.559E+06 

7.186E-02 

5.553E-02 

0.910 

52 

Airplane Mode 

70_405_TL.rpr 

0 

0.70 

40.5 

Y 

N 

1.471E+04 

2.073E+06 

9.627E-02 

7.383E-02 

0.917 

80 

Airplane Mode 

70_37_TL.rpr 

0 

0.70 

37.0 

Y 

N 

1.369E+03 

3.036E+05 

6.184E-03 

6.200E-03 

0.583 

53 

Hover 

6_TL_SD_GC.rpr 

90 

0.00 

6.0 

Y 

Y 

NA 

NA 

NA 

NA 

NA 

54 

Hover 

6_TL_SD_GE.rpr 

90 

0.00 

6.0 

Y 

Y 

NA 

NA 

NA 

NA 

NA 

55 

Hover 

6_TL_SD_IR.rpr 

90 

0.00 

6.0 

Y 

Y 

NA 

NA 

NA 

NA 

NA 

56 

Hover 

10_TL_SD_IR2.rpr 

90 

0.00 

6.0 

Y 

Y 

NA 

NA 

NA 

NA 

NA 

57 

Advance Ratio 

27_7 5_2_TL_S D . r p r 

75 

0.27 

2.0 

Y 

N 

1.812E+04 

3.915E+05 

8.185E-02 

7.995E-03 

NA 

58 

Advance Ratio 

27_7 5_3_TL_S D . r p r 

75 

0.27 

3.0 

Y 

N 

2.173E+04 

4.780E+05 

9.815E-02 

9.761E-03 

NA 

59 

Advance Ratio 

2 7_7 5_4_T L_S D . r p r 

75 

0.27 

4.0 

Y 

N 

2.530E+04 

5.765E+05 

1.143E-01 

1.177E-02 

NA 

88 

Advance Ratio 

27_7 5_1_TL_S D . r p r 

75 

0.27 

1.0 

Y 

N 

NA 

NA 

NA 

NA 

NA 

87 

Advance Ratio 

2 7_7 5_0_T L_S D . r p r 

75 

0.27 

0.0 

Y 

N 

1.089E+04 

2.507E+05 

4.919E-02 

5.120E-03 

NA 

61 

Advance Ratio 

18_7 5_3_TL_S D . r p r 

75 

0.18 

3.0 

Y 

N 

2.368E+04 

4.659E+05 

1.070E-01 

9.514E-03 

NA 

62 

Advance Ratio 

18_75_4_TL_SD.rpr 

75 

0.18 

4.0 

Y 

N 

2.681E+04 

5.488E+05 

1.211E-01 

1.121E-02 

NA 

63 

Advance Ratio 

18_75_5_TL_SD.rpr 

75 

0.18 

5.0 

Y 

N 

2.988E+04 

6.429E+05 

1.350E-01 

1.313E-02 

NA 

83 

Advance Ratio 

18_75_l_TL_SD.rpr 

75 

0.18 

1.0 

Y 

N 

1.733E+04 

3.264E+05 

7.828E-02 

6.666E-03 

NA 

84 

Advance Ratio 

18_7 5_2_TL_S D . r p r 

75 

0.18 

2.0 

Y 

N 

2.049E+04 

3.915E+05 

9.255E-02 

7.995E-03 

NA 

85 

Advance Ratio 

18_7 5_3_TL_S D . r p r 

75 

0.18 

3.0 

Y 

N 

2.368E+04 

4.659E+05 

1.070E-01 

9.514E-03 

NA 

86 

Advance Ratio 

1 8_7 5_0_T L_S D . r p r 

75 

0.18 

0.0 

Y 

N 

1.415E+04 

2.704E+05 

6.392E-02 

5.522E-03 

NA 

64 

Hover 

10_SD.rpr 

90 

0.00 

10.0 

N 

Y 

4.447E+04 

1.291E+06 

1.932E-01 

2.486E-02 

0.688 

65 

Hover 

0_SD.rpr 

90 

0.00 

0.0 

N 

Y 

1.682E+04 

2.903E+05 

7.307E-02 

5.591E-03 

0.711 

66 

Hover 

2_SD.rpr 

90 

0.00 

2.0 

N 

Y 

2.237E+04 

4.075E+05 

9.718E-02 

7.848E-03 

0.773 

67 

Hover 

4_SD.rpr 

90 

0.00 

4.0 

N 

Y 

2.800E+04 

5.617E+05 

1.216E-01 

1.082E-02 

0.790 

68 

Hover 

6_SD.rpr 

90 

0.00 

6.0 

N 

Y 

3.391E+04 

7.559E+05 

1.473E-01 

1.456E-02 

0.781 

69 

Hover 

8_SD.rpr 

90 

0.00 

8.0 

N 

Y 

3.962E+04 

9.945E+05 

1.721E-01 

1.915E-02 

0.749 

89 

Tilt 

15_18.rpr 

15 

0.32 

18.0 

N 

N 

8.523E+03 

6.435E+05 

3.850E-02 

1.314E-02 

NA 

90 

Tilt 

75_3.rpr 

75 

0.32 

3.0 

N 

N 

1.705E+04 

4.130E+05 

7.701E-02 

8.434E-03 

NA 

91 

Advance Ratio 

18_75_3.rpr 

75 

0.18 

3.0 

N 

N 

2.075E+04 

3.998E+05 

9.373E-02 

8.164E-03 

NA 

92 

Airplane Mode 

40_24.rpr 

0 

0.40 

24.0 

N 

N 

6.246E+03 

5.274E+05 

4.088E-02 

1.878E-02 

0.874 
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Appendix D Unsteady XV-1 5 Rotor Validation Results 

D. 1 Simulations Parameters 

Table D.5 shows a compact overview of the simulation parameters. ND stands for non-dimensionalized, the boundary size is 
expressed as the coordinates of the corners of a rectangular prism. The grid cells indicate the amount of cells on the x,y and 
z edge of the boundary, respectively. Rotor Grid Ref. and Grid Refinement indicate refinement of cells by multiplying the 
amount of cells by a certain factor. 

Table D.5 Overview of steady simulation parameters. 


Hover Parameters 


Boundary Conditions & 

Grid 


Simulation Variables 


Radius 

3.81 m 

Boundary Size 


(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.1) 

Simulated time 

2.5 s 

Cone angle 

2.5 deg 

Grid cells [#] 


30,30,42 

Timesteps 

3000 

Tip speed 

225.55 m/s 

Rotor Box Ref. 


1,1 

AT 

8.33E-04 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 


6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement Box [m] 


(-7.62, -7.62, -68.58) 
(7.62, 7.62, 38.1) 

Iterations 

10 

Power ND 

6.41E+08 

Grid Refinement 


4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

2.84E+06 

Volume Ratio Max 


8 

Flight Condition 

hover 



Cells 


1,053,948 

Rotor Model 


Tilt and Advance Ratio 

Boundary Conditions & 

Grid 


Simulation Variables 


Radius 

3.81 m 

Boundary Size 


(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.10) 

Simulated time 

1.25 s 

Cone angle 

2.5 deg 

Grid cells [#] 


30,30,42 

Timesteps 

1000 

Tip speed 

221.17 m/s 

Rotor Box Ref. 


1,1 

AT 

8.33E-04 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 


6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement Box [m] 


(-15.24, -7.62, -68,58) 
(7.62,7.62 38.10) 

Iterations 

10 

Power ND 

6.05E+08 

Grid Refinement 


4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

2.73E+06 

Volume Ratio Max 


8 

Flight Condition 

general 



Cells 


2,029,692 (at 30 degree tilt) 



Airplane Mode 


Boundary Conditions & 

Grid 


Simulation Variables 


Radius 

3.81 m 

Boundary Size 


(-38.1, -38.1,-68.58) 
(38.1, 38.1, 38.10) 

Simulated time 

1.25 s 

Cone angle 

2.5 deg 

Grid cells [#] 


30,30,42 

Timesteps 

1000 

Tip speed 

183.76 m/s 

Rotor Box Ref. 


1,1 

AT 

8.33E-04 

Cutout 

0.0875 r/R 

Rotor Grid Ref. 


6x 

AX min 

0.0784 

Solidity 

0.081 

Grid Refinement Box [m] 


(-15.24, -7.62, -68,58) 
(7.62,7.62 38.10) 

Iterations 

10 

Power ND 

3.47E+08 

Grid Refinement 


4x 

Relaxation (u,v,w,p) 

0.1 

Thrust ND 

1.89E+06 

Volume Ratio Max 


8 

Flight Condition 

general 


Cells 1,053,948 
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D.2 Hover - Unsteady 

Figure D.10 and Figure D.1 1 show, analogous to the steady results obtained in Section 4.1, the hover performance for the 
unsteady model of the XV-15 rotor in RotUNS. 
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Figure D.1 0 Unsteady results for XV- 1 5 rotor hover power as a function of thrust [2], [28]. 
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Both graphs show serious error in excess of 1 0%, compared to the presented theoretical and experimental data. 


1.000 


0.800 


0.600 


0 

E 

o 

2 

3 

q= 


0.400 


0.200 


0.000 


0.00 





4 



♦ 


▲ 

♦ 


OARF Data (Mtip = 0.66) 
WTT (Mtip = .53) 


CAMRAD I (Mtip = .66) 

RotCFD UNSTEADY TL (Mtip = .66) 

O RotCFD UNSTEADY TL SD E2 (Mtip = .66) 
O RotCFD UNSTEADY TL SD E4 (Mtip = .66) 
CAMRAD II (Mtip = .69) 


I Q n 





WTT (Mtip = .65) 

Flight test (Mtip = .69) 

RotCFD UNSTEADY (Mtip = .66) 

RotCFD UNSTEADY TL SD (Mtip = .66) 

□ RotCFD UNSTEADY TL SD E3 (Mtip = .66) 
RotCFD UNSTEADY SD (Mtip = .66) 



0.05 0.10 0.15 0.20 

blade loading, C T /o [~] 

Figure D.1 1 Unsteady results for XV-15 rotor hover figure of Merit as a function of thrust [2], [28]. 


D.3 Tilt Mode - Unsteady 

Figure D.1 2 and Figure D.1 3 show, analogous to the steady results obtained in Section 4.2, the tilt mode performance and 
the sensitivity to advance ratio variations for the unsteady model of the XV-1 5 rotor in RotUNS. 
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Figure D.12 Unsteady results for XV-15 rotor power as a function of thrust for various pylon angles at V/QR = .32 [2]. 

The results for the various tilt modes correlate slightly better than the hover performance, the steady model, however, 
outperforms the unsteady model in terms of accuracy. 
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Figure D.13 Unsteady results for XV-15 rotor power as function of thrust, for a P = 75° and Mtip = 0.65 [2]. 
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The variation of advance ratio shows similar performance to CAMRAD I. 
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D.4 Airplane Mode - Unsteady 

Figure D.14 shows, analogous to the steady results obtained in Section 4.3, airplane mode performance for the unsteady 
model of the XV-15 rotor in RotUNS. 
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Figure D.14 Unsteady results for (rotor) propulsive efficiency as function of thrust [2]. 
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D.5 Data Unsteady XV- 1 5 Rotor Results 

The data used in the various plots is tabulated in Table D.6. 


Table D.6 Summary of unsteady RotCFD Validation Data. 


# 

Flight Mode 

File 

a P 

[deg] 

V/OR 

[m/s] 

00 

[deg] 

TL 

SD 

T [N] 

P [J/s] 

C t/o [-] 

Cp/G [-] 

M [-] or 

nH 

1 

Hover 

O.rpr 

90 

0.00 

0.0 

N 

N 

1.100E+04 

2.170E+05 

4.778E-02 

4.179E-03 

0.505 

2 

Hover 

2.rpr 

90 

0.00 

2.0 

N 

N 

1.460E+04 

3.020E+05 

6.342E-02 

5.816E-03 

0.552 

3 

Hover 

4.rpr 

90 

0.00 

4.0 

N 

N 

1.820E+04 

4.130E+05 

7.906E-02 

7.954E-03 

0.562 

4 

Hover 

6.rpr 

90 

0.00 

6.0 

N 

N 

2.170E+04 

5.510E+05 

9.427E-02 

1.061E-02 

0.552 

5 

Hover 

8.rpr 

90 

0.00 

8.0 

N 

N 

2.500E+04 

7.170E+05 

1.086E-01 

1.381E-02 

0.524 

6 

Hover 

lO.rpr 

90 

0.00 

10.0 

N 

N 

2.800E+04 

9.080E+05 

1.216E-01 

1.749E-02 

0.489 

7 

Hover 

0_TL.rpr 

90 

0.00 

0.0 

Y 

N 

1.070E+04 

2.160E+05 

4.648E-02 

4.160E-03 

0.481 

8 

Hover 

2_TL.rpr 

90 

0.00 

2.0 

Y 

N 

1.400E+04 

3.000E+05 

6.082E-02 

5.778E-03 

0.526 

9 

Hover 

4_TL.rpr 

90 

0.00 

4.0 

Y 

N 

1.730E+04 

4.070E+05 

7.515E-02 

7.839E-03 

0.531 

10 

Hover 

6_TL.rpr 

90 

0.00 

6.0 

Y 

N 

2.040E+04 

5.410E+05 

8.862E-02 

1.042E-02 

0.513 

11 

Hover 

8_TL.rpr 

90 

0.00 

8.0 

Y 

N 

2.350E+04 

7.100E+05 

1.021E-01 

1.367E-02 

0.483 

12 

Hover 

10_TL.rpr 

90 

0.00 

10.0 

Y 

N 

2.620E+04 

9.070E+05 

1.138E-01 

1.747E-02 

0.442 

13 

Hover 

0_TL_SD.rpr 

90 

0.00 

0.0 

Y 

Y 

1.270E+04 

2.540E+05 

5.517E-02 

4.892E-03 

0.545 

14 

Hover 

2_TL_SD.rpr 

90 

0.00 

2.0 

Y 

Y 

1.698E+04 

3.492E+05 

7.375E-02 

6.726E-03 

0.584 

15 

Hover 

4_TL_SD.rpr 

90 

0.00 

4.0 

Y 

Y 

2.020E+04 

4.750E+05 

8.775E-02 

9.148E-03 

0.573 

16 

Hover 

6_TL_SD.rpr 

90 

0.00 

6.0 

Y 

Y 

2.390E+04 

6.310E+05 

1.038E-01 

1.215E-02 

0.554 

17 

Hover 

8_TL_SD.rpr 

90 

0.00 

8.0 

Y 

Y 

2.730E+04 

8.170E+05 

1.186E-01 

1.574E-02 

0.520 

18 

Hover 

10_TL_SD.rpr 

90 

0.00 

10.0 

Y 

Y 

3.040E+04 

1.040E+06 

1.321E-01 

2.003E-02 

0.480 
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19 

Tilt 

15_18_TL.rpr 

15 

0.32 

18.0 

Y 

N 

4.660E+03 

3.970E+05 

2.105E-02 

8.107E-03 

NA 

20 

Tilt 

15_19_TL.rpr 

15 

0.32 

19.0 

Y 

N 

7.020E+03 

5.770E+05 

3.171E-02 

1.178E-02 

NA 

21 

Tilt 

15_20_TL.rpr 

15 

0.32 

20.0 

Y 

N 

9.390E+03 

7.670E+05 

4.241E-02 

1.566E-02 

NA 

22 

Tilt 

15_21_TL.rpr 

15 

0.32 

21.0 

Y 

N 

1.160E+04 

9.610E+05 

5.240E-02 

1.962E-02 

NA 

23 

Tilt 

30_155_TL.rpr 

30 

0.32 

15.5 

Y 

N 

4.918E+03 

3.443E+05 

2.221E-02 

7.031E-03 

NA 

24 

Tilt 

30_165_TL.rpr 

30 

0.32 

16.5 

Y 

N 

7.218E+03 

5.011E+05 

3.260E-02 

1.023E-02 

NA 

25 

Tilt 

30_175_TL.rpr 

30 

0.32 

17.5 

Y 

N 

9.489E+03 

6.648E+05 

4.286E-02 

1.358E-02 

NA 

26 

Tilt 

30_185_TL.rpr 

30 

0.32 

18.5 

Y 

N 

1.175E+04 

8.360E+05 

5.307E-02 

1.707E-02 

NA 

27 

Tilt 

60_8_TL.rpr 

60 

0.32 

8.0 

Y 

N 

8.855E+03 

3.655E+05 

4.000E-02 

7.463E-03 

NA 

28 

Tilt 

60_9_TL.rpr 

60 

0.32 

9.0 

Y 

N 

1.104E+04 

4.647E+05 

4.985E-02 

9.489E-03 

NA 

29 

Tilt 

60_10_TL.rpr 

60 

0.32 

10.0 

Y 

N 

1.316E+04 

5.738E+05 

5.944E-02 

1.172E-02 

NA 

30 

Tilt 

60_ll_TL.rpr 

60 

0.32 

11.0 

Y 

N 

1.524E+04 

6.855E+05 

6.884E-02 

1.400E-02 

NA 

31 

Tilt 

75_2_TL.rpr 

75 

0.32 

2.0 

Y 

N 

8.970E+03 

2.660E+05 

4.052E-02 

5.432E-03 

NA 

32 

Tilt 

75_3_TL.rpr 

75 

0.32 

3.0 

Y 

N 

1.120E+04 

3.250E+05 

5.059E-02 

6.637E-03 

NA 

33 

Tilt 

75_4_TL.rpr 

75 

0.32 

4.0 

Y 

N 

1.330E+04 

3.920E+05 

6.008E-02 

8.005E-03 

NA 

34 

Tilt 

75_5_TL.rpr 

75 

0.32 

5.0 

Y 

N 

1.543E+04 

4.673E+05 

6.970E-02 

9.544E-03 

NA 

35 

Advance Ratio 

75_2_TL.rpr 

75 

0.32 

2.0 

Y 

N 

8.970E+03 

2.660E+05 

4.052E-02 

5.432E-03 

NA 

36 

Advance Ratio 

75_3_TL.rpr 

75 

0.32 

3.0 

Y 

N 

1.120E+04 

3.250E+05 

5.059E-02 

6.637E-03 

NA 

37 

Advance Ratio 

75_4_TL.rpr 

75 

0.32 

4.0 

Y 

N 

1.330E+04 

3.920E+05 

6.008E-02 

8.005E-03 

NA 

38 

Advance Ratio 

75_5_TL.rpr 

75 

0.32 

5.0 

Y 

N 

1.543E+04 

4.673E+05 

6.970E-02 

9.544E-03 

NA 

39 

Advance Ratio 

27_75_2_TL.rpr 

75 

0.27 

2.0 

Y 

N 

1.007E+04 

2.710E+05 

4.548E-02 

5.533E-03 

NA 

40 

Advance Ratio 

27_75_3_TL.rpr 

75 

0.27 

3.0 

Y 

N 

1.220E+04 

3.260E+05 

5.511E-02 

6.657E-03 

NA 

41 

Advance Ratio 

27_75_4_TL.rpr 

75 

0.27 

4.0 

Y 

N 

1.424E+04 

3.911E+05 

6.431E-02 

7.987E-03 

NA 

42 

Advance Ratio 

18_75_3_TL.rpr 

75 

0.18 

3.0 

Y 

N 

1.390E+04 

3.330E+05 

6.279E-02 

6.800E-03 

NA 

43 

Advance Ratio 

18_75_4_TL.rpr 

75 

0.18 

4.0 

Y 

N 

1.577E+04 

3.913E+05 

7.123E-02 

7.991E-03 

NA 

44 

Advance Ratio 

18_75_5_TL.rpr 

75 

0.18 

5.0 

Y 

N 

1.760E+04 

4.560E+05 

7.950E-02 

9.312E-03 

NA 

45 

Airplane Mode 

40_23_TL.rpr 

0 

0.40 

23.0 

Y 

N 

1.404E+03 

1.670E+05 

6.341E-03 

3.410E-03 

0.621 

46 

Airplane Mode 

40_24_TL.rpr 

0 

0.40 

24.0 

Y 

N 

3.128E+03 

3.016E+05 

1.413E-02 

6.159E-03 

0.766 

47 

Airplane Mode 

40_25_TL.rpr 

0 

0.40 

25.0 

Y 

N 

4.820E+03 

4.385E+05 

2.177E-02 

8.955E-03 

0.811 

48 

Airplane Mode 

40_26_TL.rpr 

0 

0.40 

26.0 

Y 

N 

6.478E+03 

5.794E+05 

2.926E-02 

1.183E-02 

0.825 

49 

Airplane Mode 

70_375_TL.rpr 

0 

0.70 

37.5 

Y 

N 

1.397E+03 

3.099E+05 

6.308E-03 

6.329E-03 

0.582 

50 

Airplane Mode 

70_385_TL.rpr 

0 

0.70 

38.5 

Y 

N 

3.607E+03 

6.043E+05 

1.629E-02 

1.234E-02 

0.771 

51 

Airplane Mode 

70_395_TL.rpr 

0 

0.70 

39.5 

Y 

N 

5.737E+03 

9.014E+05 

2.591E-02 

1.841E-02 

0.822 

52 

Airplane Mode 

70_405_TL.rpr 

0 

0.70 

40.5 

Y 

N 

7.800E+03 

1.200E+06 

3.523E-02 

2.451E-02 

0.840 

53 

Hover 

6_TL_SD_GC.rpr 

90 

0.00 

6.0 

Y 

Y 

2.390E+04 

6.159E+05 

1.038E-01 

1.186E-02 


54 

Hover 

6_TL_SD_GE.rpr 

90 

0.00 

6.0 

Y 

Y 

2.390E+04 

6.300E+05 

1.038E-01 

1.213E-02 

0.552 

55 

Hover 

6_TL_SD_IR.rpr 

90 

0.00 

6.0 

Y 

Y 

2.359E+04 

6.299E+05 

1.025E-01 

1.213E-02 

0.545 

56 

Hover 

10_TL_SD_IR2.rpr 

90 

0.00 

6.0 

Y 

Y 

2.980E+04 

1.030E+06 

1.295E-01 

1.984E-02 

0.472 

58 

Advance Ratio 

27_7 5_2_TL_S D . r p r 

75 

0.27 

2.0 

Y 

N 

1.202E+04 

3.146E+05 

5.428E-02 

6.425E-03 

NA 

59 

Advance Ratio 

27_7 5_3_TL_S D . r p r 

75 

0.27 

3.0 

Y 

N 

1.435E+04 

3.789E+05 

6.482E-02 

7.738E-03 

NA 

60 

Advance Ratio 

27_75_4_TL_SD.rpr 

75 

0.27 

4.0 

Y 

N 

1.669E+04 

4.516E+05 

7.539E-02 

9.222E-03 

NA 

61 

Advance Ratio 

18_75_3_TL_SD.rpr 

75 

0.18 

3.0 

Y 

N 

1.619E+04 

3.868E+05 

7.313E-02 

7.899E-03 

NA 

62 

Advance Ratio 

1 8_7 5 _4_T L_S D . r p r 

75 

0.18 

4.0 

Y 

N 

1.832E+04 

4.538E+05 

8.275E-02 

9.267E-03 

NA 

63 

Advance Ratio 

18_7 5_5_TL_S D . r p r 

75 

0.18 

5.0 

Y 

N 

2.037E+04 

5.276E+05 

9.201E-02 

1.077E-02 

NA 
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64 

Hover 

10_SD.rpr 

90 

0.00 

10.0 

N 

Y 

3.267E+04 

1.044E+06 

1.419E-01 

2.010E-02 

0.526 

65 

Hover 

0_SD.rpr 

90 

0.00 

0.0 

N 

Y 

1.278E+04 

2.525E+05 

5.552E-02 

4.863E-03 

0.526 

66 

Hover 

2_SD.rpr 

90 

0.00 

2.0 

N 

Y 

1.696E+04 

3.515E+05 

7.368E-02 

6.769E-03 

0.608 

67 

Hover 

4_SD.rpr 

90 

0.00 

4.0 

N 

Y 

2.118E+04 

4.801E+05 

9.201E-02 

9.246E-03 

0.601 

68 

Hover 

6_SD.rpr 

90 

0.00 

6.0 

N 

Y 

2.512E+04 

6.389E+05 

1.091E-01 

1.231E-02 

0.593 

69 

Hover 

8_SD.rpr 

90 

0.00 

8.0 

N 

Y 

2.903E+04 

8.273E+05 

1.261E-01 

1.593E-02 

0.563 
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Appendix E Wind Tunnel Case Plots 

Each of the wind tunnel subsets would have at least a couple of interesting velocity or vector plots to show, however, since 
a total of 28 subsets have been processed this would pose too much images in this Appendix. Therefore only a selection of 
view per subset is presented here. The used coordinate system in the NFAC wind tunnels is shown in Figure E.15. The TTR 
will always rotate for the thrust to be (partially) aligned with the positive y-axis. 




Figure E.15 The boundaries of the extended test sections with TTR in edgewise and axial mode, respectively. 

E. 1 80- by 120-Foot Wind Tunnel Cases 

E.1.1 Case 1, Edgewise, 0 kts. 



Figure E.16 Case 1, WTGE, ZY-plane. Figure E.17 Case 1, WTRO, ZY-plane. 

Figure E.16 and Figure E.17 show a velocity plot (with pressure plot on the geometry) of the ZY-plane of case 1 for WTGE 
and WTRO, respectively. Figure E.1 8 serves as the legend for the plots of case 1 , all units are SI. 



Pressure, p [Pa] : 1.020E+05 1.022E+05 1.024E+05 1.026E+05 1.028E+05 1.031 E+05 1.033E+05 1.035E+05 


Velocity Magnitude, V [m/s]: 0.00E+00 4.08E+00 8.16E+00 1.22E+01 1.63E+01 2.04E+01 2.45E+01 2.86E+01 3.27E+01 3.67E+01 4.08E+01 4.49E+01 4.90E+01 

Figure E.1 8 Legend for pressure and velocities for case 1 . 
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E. 1.2 Case 2, Edgewise, lOkts. 



Figure E.19 Case 2, FFGE, XY-plane. Figure E.20 Case 2, WTGE, XY-plane. 

Figure E.19 and Figure E.20 show a velocity plot (with pressure plot on the geometry) of the XY-plane of case 2 for FFGE and 
WTGE, respectively. Figure E.21 serves as the legend for the plots of case 2. 



Pressure, p [Pa] : 1 .020E+05 1 .022E+05 1 .024E+05 1 .026E+05 


Velocity Magnitude, V [m/s]: 0.00E+00 4.08E+00 8.16E+00 1.22E+01 1.63E+01 2.04E+01 2.45E+01 2.86E+01 3.27E+01 3.67E+01 4.08E+01 4.49E+01 

Figure E.21 Legend for pressure and velocities for case 2. 


1 .035E+05 
4.90E+01 
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E. 1.3 Case 3, Axial, 1 00 kts. 



Figure E.22 Case 3, WTGE, XZ-plane. Figure E.23 Case 3, WTRO, XZ-plane. 


Figure E.22 and Figure E.23 show a velocity plot (with pressure plot on the geometry) of the ZY-plane of case 3 for WTGE 
and WTRO, respectively. Figure E.24 serves as the legend for the plots of case 3. 


p [Pa] : 9.69E+04 9.73E+04 9.76E+04 9.80E+04 9.84E+04 9.87E+04 9.91E+04 9.95E+04 9.98E+04 1.00E+05 1.01E+05 1.01E+05 1.01E+05 1.02E+05 1.02E+05 1.02E+05 1.03E+05 
V [m/s]: 0.00E+00 5.11E+00 1.02E+01 1.53E+01 2.05E+01 2.56E+01 3.07E+01 3.58E+01 4.09E+01 4.60E+01 5.11E+01 5.62E+01 6.14E+01 6.65E+01 7.16E+01 7.67E+01 8.18E+01 

Figure E.24 Legend for pressure and velocities for case 3. 
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E.2 40- by 80-Foot Wind Tunnel Cases 

E.2 . 7 Case 5, Edgewise, 1 00 kts. 



Figure E.25 Case 5, FFGE, XY-plane. Figure E.26 Case 5, WTGE, XY-plane. 

Figure E.25 and Figure E.26 show a velocity plot (with pressure plot on the geometry) of the XY-plane of case 5 for FFGE and 
WTGE, respectively. Figure E.27 serves as the legend for the plots of case 5. 



Pressure, p [Pa]: 9.674E+04 9.748E+04 9.821E+04 9.895E+04 9.969E+04 1.004E+05 1.012E+05 1.019E+05 1.026E+05 

Velocity Magnitude, V [m/s]: 0.00E+00 7.63E+00 1.53E+01 2.29E+01 3.05E+01 3.82E+01 4.58E+01 5.34E+01 6.10E+01 6.87E+01 7.63E+01 8.39E+01 9.16E+01 9.92E+01 

Figure E.27 Legend for pressure and velocities for case 5. 
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E.2.2 Case 6, Axial, 200 kts. 



Figure E.28 Case 6, FFGE, XY-plane. Figure E.29 Case 6, WTGE, XY-plane. 

Figure E.28 and Figure E.29 show a velocity plot (with pressure plot on the geometry) of the XY-plane of case 6 for FFGE and 
WTGE, respectively. Figure E.30 serves as the legend for the plots of case 6. 


Pressure, p [Pa] : 9.000E+04 9.200E+04 9.400E+04 9.600E+04 9.800E+04 1.000E+05 1.020E+05 1.040E+05 

Velocity Magnitude, V [m/s]: 0.00E+00 1.41E+01 2.82E+01 4.23E+01 5.65E+01 7.06E+01 8.47E+01 9.88E+01 1.13E+02 1.27E+02 1.41E+02 1.55E+02 1.69E+02 1.84E+02 

Figure E.30 Legend for pressure and velocities for case 6. 

E.2.3 Case 7 , Tilt, 100 kts. 



Figure E.31 Case 7, WTGE, XY-plane. Figure E.32 Case 7, WTRO, XY-plane. 


Figure E.31 and Figure E.32 show a velocity plot (with pressure plot on the geometry) of the XY-plane of case 7 for WTGE 
and WTRO, respectively. Figure E.33 serves as the legend for the plots of case 7. The grid at the walls of the 40-by 80-Foot 
Wind Tunnel shows problems with gridding clearly, as shown in Figure E.34. 


Pressure, p [Pa] : 8.000E+04 8.353E+04 


8.706E+04 9.059E+04 9.412E+04 9.765E+04 1.012E+05 1.047E+05 1.082E+05 


Velocity Magnitude, V [m/s]: 0.00E+00 8.03E+00 1.61E+01 2.41E+01 3.21E+01 4.02E+01 4.82E+01 5.62E+01 6.43E+01 7.23E+01 8.03E+01 8.84E+01 9.64E+01 

Figure E.33 Legend for pressure and velocities for case 7. 


Wind Tunnel Interference Effects on Tilt Rotor Testing Using Computational Fluid Dynamics Witold J. F. Koning 


81 





Figure E.34 ZY-plane of case 7, just behind the rotor plane. 
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